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SECTION  I 
INTRODUCTION 

The  numerical  methods  for  the  computation  of  aerodynamic  forces  of  small- 
disturbance  transonic  flows  about  oscillating  airfoils  and  planar  wings  have 
been  rapidly  developed  in  recent  years.  Such  developments  (References  1-32) 
are  briefly  discussed  In  the  next  section.  The  numerical  methods  for  the 
computation  of  free  vibration  frequencies  and  normal  modes  and  the  time- 
history  dynamic  responses  of  wing  structures  have  been  well -developed.  It 
Is  the  time  to  integrate  the  developments  in  transonic  aerodynamics  and 
structural  dynamics  for  aeroelastic  applications. 

Such  aeroelastic  applications  have  recently  begun.  In  References  5, 
Traci,  et_  al_.  performed  a flutter  analysis  for  the  NACA  64A006  and  NACA 
64A410  airfoils  based  on  the  aerodynamic  coefficients  computed  by  using 
relaxation  programs  STRANS  and  UTRANS  (References  4 and  6).  The  aerodynamic 
coefficients  were  obtained  for  the  case  that  the  airfoil  pitches  about  the 
leading  edge.  The  moment  coefficients  were  then  transformed  from  the 
leading  edge  to  the  elastic  axis  for  flutter  analysis.  The  flutter  velocities 
and  frequencies  were  obtained  by  using  the  U-g  method  (Reference  33).  Primary 
results  were  obtained  for  the  NACA  64A006  airfoil  with  two  degrees  of 
freedom,  plunge  and  pitch.  The  results  demonstrated  the  crucial  importance 
of  the  transonic  range  by  way  of  the  drop  in  the  flutter  velocity  parameter 
at  M = 0.85  and  the  corresponding  large  variation  in  flutter  frequency  in 
this  range.  The  restricted  Mach  number  range  of  the  aerodynamic  data  did  not 
permit  an  in-depth  examination  of  NACA  64A410  airfoil. 

The  effect  of  mean  angle  of  attack  and  wind  tunnel  wall  interference  on 
the  flutter  parameters  for  the  NACA  64A006  airfoil  was  also  examined  for  the 
limited  case  for  which  the  aerodynamic  coefficients  had  been  calculated. 
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Flutter  calculations  were  also  attempted  for  the  NACA  64A006  airfoil 
with  an  additional  degree  of  freedom  of  aileron  oscillation.  Due  to  the 
rather  restricted  reduced  frequency  range  of  their  aerodynamic  data,  the  true 
flutter  velocity  was  not  obtainable  because  an  aileron  dominated  mode  was 
fluttering  at  a frequency  out  of  the  data  range.  It  was  determined,  however, 
that  the  aileron  oscillation  significantly  damped  the  plunge  mode  which 
dominated  the  flutter  boundary  in  the  two  degree  of  freedom  case.  The  curious 
coupling  of  the  aileron  oscillation  and  plunge  modes  and  the  high  frequency 
aileron  buzz  phenomena  appear  to  be  peculiar.  A thorough  examination  of  the  three 
degree  of  freedom  case,  however,  requires  the  generation  of  unsteady  aerodynamic 
coefficients  over  a wider  range  of  reduced  frequencies  than  those  calculated 
in  Reference  5. 

As  an  alternative  but  more  complete  solution,  Ballhaus  and  Goorjian 
(Reference  28)  proposed  to  use  the  time-integration  and  the  indicial  approaches 
for  the  aeroelastic  response  problems.  They  performed  an  aeroelastic  response 
analysis  for  an  NACA  64A006  airfoil  with  a single  pitching  degree  of  freedom. 

The  time-history  responses  were  computed  using  the  computer  program  LTRAN2 
for  unsteady  transonic  flow  coupled  with  an  integration  procedure  for  the 
simple  differential  equation  of  motion  for  the  airfoil.  The  motion  of  the 
airfoil  in  a flow  at  = 0.88  was  forced  for  the  first  few  cycles  with  an 
amplitude  of  0.5  degrees  until  the  pitching  moment  became  periodic,  after 
which  the  airfoil  motion  and  the  aerodynamic  response  were  left  free  to  drive 
each  other.  As  the  structural  damping  coefficient  was  varied  parametrically, 
the  responses  were  obtained  for  the  highly  damped,  neutral,  slightly  diver- 
gent, and  highly  divergent  cases.  Another  case  with  an  initial  amplitude  of 
1.5°  and  * 0.87  was  also  calculated.  This  example  illustrated  the  non- 
sinusoidal  pitching-moment  behavior  that  could  result  from  the  large  shock- 


wave  excursions  encountered  at  larger  airfoil  motion  amplitudes. 

From  the  development  in  Reference  28,  it  appears  promising  that  the 
simultaneous  integration  of  the  flow  and  structure  may  lead  to  the 
solution  of  the  transonic  aeroelastic  response  problems  for  the  multiple 
degree  of  freedom  system. 

As  a visiting  scientist  at  the  Air  Force  Flight  Dynamics  Laboratory, 
Rizzetta  published  three  reports  (References  34-36)  on  his  investigation 
of  the  transonic  aeroelastic  problems  of  two-dimensional  airfoils. 

In  Reference  34,  Rizzetta  performed  a flutter  analysis  of  an  NACA 
64A010  airfoil  with  two  degrees  of  freedom,  plunge  and  pitch.  The  steady 
and  unsteady  transonic  flows  were  computed  by  using  the  relaxation  computer 
programs  STRANS  and  UTRANS,  respectively  (Reference  4).  The  flutter 
velocities  and  frequencies  were  computed  using  the  standard  U-g  method. 

Flutter  data  were  obtained  for  two  Mach  numbers  (0.72  and  0.80)  and  two 
angles  of  attack  (0°  and  1°). 

It  was  reported  that  for  both  angles  of  attack  and  both  Mach  numbers 
the  flutter  velocity  parameter  decreased  monotonically  with  increasing 
altitude  and  with  increasing  first  mass  moment.  The  flutter  velocity 
parameter  increased  with  increasing  frequency  ratio,  defined  as  the  square 
of  the  ratio  of  the  plunge  frequency  to  the  pitch  frequency.  The  increase 
in  Mach  number  noticeably  decreased  the  flutter  speed  parameter,  especially 
in  the  presence  of  a shock.  The  study  of  the  first  mass  moment  was  for  its 
effect  on  the  flutter  suppression. 

In  Reference  35,  Rizzetta  performed  a comparative  study  of  the  relaxation 
programs  STRANS2  and  UTRANS2  and  the  time- integration  program  LTRAN2  for 


calculating  steady  and  unsteady  aerodynamic  pressures  for  an  NACA  64A010  airfoil 
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A finite  difference  mesh  of  43x70  was  used  In  the  relaxation  method  and  a 
mesh  of  79x99  was  used  in  the  time  integration  method.  Mach  numbers  considered 
were  0.72  and  0.82.  The  airfoil  was  pitching  at  an  amplitude  of  a=l°  with 
reduced  frequency  kc  of  0.05  and  0.2,  respectively.  A mean  angle  of  attack  of 
1 degree  was  also  considered  for  M = 0.72,  a = 0.5°  and  kc  = 0.05.  The 

agreement  between  the  two  solutions  was  good  except  for  the  unsteady  results 

at  M = 0.82  with  the  presence  of  shock.  In  that  case,  the  peak  lift  value 
obtained  from  the  time-integration  method  was  approximately  50%  higher  than 
that  from  the  harmonic  analysis. 

In  Reference  36,  Rizzetta  performed  an  aeroelastic  response  study  of  an 
NACA  64A010  airfoil  by  simultaneously  integrating  the  LTRAN2  aerodynamics 
program  and  the  structural  equations  of  motion.  The  study  included  a single 
degree  of  freedom  case  with  airfoil  pitching  at  quarter  chord  and  a three 

degree  of  freedom  case  with  airfoil  pitching,  plunging,  and  aileron  pitching. 

The  study,  however,  included  neither  the  single  degree  of  freedom  case  of 
airfoil  plunging  nor  the  two  degree  of  freedom  case  with  airfoil  pitching 
and  plunging. 

Mach  numbers  considered  were  0.72  and  0.82.  Results  were  obtained  as 
plots  of  nondimensional  pitching  displacement  (and  moment),  plunging 
displacement  (and  lift),  and  aileron  displacement  (and  moment)  versus  time 
for  various  values  of  the  airfoil-air  mass  ratio  y.  All  these  displacements 
and  force  quantities  in  the  figures  were  divided  by  a'(o),  the  nondimensional 
initial  pitching  velocity  defined  in  Reference  36.  The  values  used  for  a'(o) 
ranges  from  1 to  a very  high  value  of  7.5.  By  varying  the  values  of  y, 
convergent,  divergent,  and  neutral  stable  response  curves  were  obtained. 

In  this  report,  flutter  analyses  were  performed  for  an  NACA  64A006  and 
an  NACA  64A010  airfoil  with  pitching  and  plunging  degrees  of  freedom.  Two 
parallel  sets  of  aerodynamic  coefficients  were  obtained  by  using  both  the 
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relaxation  method  (STRANS2  and  UTRANS2)  and  the  Indlclal  method  (LTRAN2). 

The  two  sets  of  aerodynamic  coefficients  and  the  two  corresponding  sets  of 
flutter  results  are  compared  and  discussed. 

To  evaluate  the  linear  solution  obtained  by  using  LTRAN2  and  UTRANS2,  a 
flat  plate  pitching  about  the  quarter  chord  axis  and  plunging  In  the  flow 
with  M = 0.7  was  considered.  The  solutions  obtained  for  the  unsteady  pressure 
distributions  by  both  methods  were  found  to  be  in  close  agreement  with  the 
linear  theory  solution  based  on  Kernel  function  method  (The  computer  program 
was  provided  by  Dr.  Samuel  R.  Bland  of  NASA-Langley  through  private  communication) 
Flutter  results  corresponding  to  these  three  sets  of  aerodynamic  coefficients 
were  obtained  and  compared. 

For  the  NACA  64A006  airfoil,  the  Mach  numbers  considered  were  0.7,  0.8, 
and  0.85,  respectively.  For  the  NACA  64A010  airfoil,  the  Mach  numbers 
considered  were  0.72,  0.76,  and  0.80,  respectively.  For  each  airfoil 
and  each  Mach  number,  flutter  velocities  and  the  corresponding  reduced 
frequencies  were  obtained  by  varying  the  airfoil-air  mass  ratio  u,  plunge- 
pitch  frequency  ratio  location  of  the  mass  center  xa>  and  location  of 

the  elastic  axis  a^. 

The  two  sets  of  flutter  velocity  curves,  one  based  on  the  aerodynamic 
coefficients  computed  by  using  LTRAN2  and  the  other  computed  by  using  UTRANS2 
are  compared  and  discussed. 
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SECTION  II 


DEVELOPMENTS  IN  TRANSONIC  FLOW  COMPUTATION 


Due  to  the  fast  development  in  digital  computers,  advances  in  the 

numerical  computations  of  transonic  flow  fields  around  oscillating  two- 

dimensional  airfoils  have  been  extensive. 

The  simplified  basic  aerodynamic  equations,  following  the  assumptions 

that  the  flow  is  two-dimensional,  inviscid,  transonic  (M  ^1),  and  that  the 

00 

velocity  disturbances  are  small  compared  to  the  free  stream  velocity 
U,  can  be  deduced  from  the  general  equation  of  continuity  of  gas  dynamics 
as  (Reference  1), 

kc2  "At  ♦ 2kc  M^xt  ■ Vc  *xx  + *yy  (1) 

where  k£  = wc/U  is  the  reduced  frequency;  M is  the  free  stream  Mach  number; 

9 m 

<f>  is  the  disturbance  velocity  potential;  V = 1 - M - (y  +1)  M <j>  ; m is  a 

c ^ 

function  of  M;  and  y is  the  ratio  of  specific  heats. 

In  deriving  the  above  equation,  the  coordinate  system  is  fixed  with 

respect  to  the  airfoil,  and  x is  aligned  with  the  free  stream  direction. 

The  flow  is  defined  as  locally  subsonic  or  supersonic,  relative  to  the 

fixed  coordinate  system,  for  V >0  or  V <0,  respectively.  A measure  of  the 

^ c 

degree  of  unsteadiness  is  given  by  the  reduced  frequency  kc  when  the  airfoil 
is  oscillating  periodically  with  a frequency  w. 

Several  numerical  approaches  have  been  developed  for  the  solution  of 
transonic  flow  fields  governed  by  Equation  1. 


(a)  The  Harmonic  or  Relaxation  Approach 


This  appraoch  was  first  developed  for  steady  flows  by  Murman  and  Cole 
(Reference  2).  It  was  improved  and  extended  by  many  computational  aerodynami cists 
to  treat  the  small  - disturbance  steady  and  unsteady  transonic  flows  about 
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oscillating  airfoils  and  planar  wings  (see,  for  example.  References  2 to  32). 

To  solve  for  Equation  1,  It  requires  the  solution  of  the  following  perturbation 
equation  for  the  potential  functions 

<fr(x,y,z,t)  * <|>0(x,y,z)  + e<h (x,y,z)e1wt  (2) 

where  $ Is  the  small-disturbance  velocity  potential;  e « 1 is  related  to  the 

small  amplitude  of  the  body  motion;  <J>o  is  the  mean  steady-state  solution 

and  <t>i  is  dependent  upon  <|>o.  The  small -disturbance  potential  equation  for  the 

steady  flow  Is  nonlinear.  The  equation  for  the  first  order  unsteady  flow  Is 
2 

linear  when  $ is  dropped,  and  for  the  case  of  harmonic  boundary  disturbances. 

Its  reduced  wave  equation  form  Is  of  mixed  elliptical /hyperbolic  type  depending 
upon  the  nature  of  the  steady  solution.  Both  steady  and  unsteady  equations  can 
be  solved  by  the  mixed  differencing  line  relaxation  procedure  first  proposed 
for  steady  flows  by  Murman  and  Cole  (Reference  2). 

Based  on  this  approach,  Traci  et  al_.  developed  four  computer  programs: 
STRANS;  UTRANS;  TDSTRN ; and  TDUTRN . Programs  STRANS  and  UTRANS  are  for  two- 
dimensional  small-disturbance  steady  and  unsteady  transonic  flows  about 
oscillating  airfoils,  respectively  (References  3 and  4).  Programs  TDSTRN 
and  TDUTRN  are  the  three-dimensional  versions  of  STRANS  and  UTRANS,  respectively, 
for  flows  about  oscillating  rectangular  planar  wings  (References  5 and  6). 

A user’s  guide  for  the  improved  version  of  the  computer  programs  STRANS2  and 
UTRANS2  was  prepared  by  Holman  (Reference  8). 

Within  the  basic  framework  of  the  relaxation  method,  Ehlers  (Reference  9) 
and  Weatherill,  Sebastian,  and  Ehlers  (Reference  10)  developed  numerical 
methods  comparable  to  those  done  In  References  3 to  8 but  less  extensive. 

They  computed  the  flow  fields  around  a flat  plate,  an  NACA  64A006  airfoil, 
and  a rectangular  planform  wing  with  an  aspect  ratio  of  5.  Ballhaus  and  Bailey 


(Reference  11)  dealt  with  a constant  chord,  23.75°  sweptback  planform 
model  with  an  aspect  ratio  of  4 and  with  a Lockheed  Cl 41  airfoil  section. 
Ballhaus,  Bailey,  and  Frick  later  included  two  additional  terms  in  the  governing 
flow  equations  and  computed  the  surface  pressures  for  the  HiMAT  wing  and 
the  ONERA  M-6  wing  (Reference  12).  Krupp  and  Cole  (Reference  13)  analyzed 
an  NACA  0012  airfoil  performing  plunging  oscillations  at  high  subsonic 
speeds.  The  results  showed  the  existence  of  a moving  shock  in  the  locally 
supersonic  zone. 

In  a recent  study  (Reference  14),  Ballhaus,  Jameson,  and  Albert  compared 
the  implicit  approximate- factorization  (AF)  algorithms  with  the  successive 
line  over-relaxation  (SLOR)  method  for  the  solution  of  nonlinear,  two- 
dimensional,  steady  transonic  small  disturbance  equation  for  flows  about 
several  airfoils.  They  concluded  that  the  AF  method  required  substantially 
less  computer  time  than  the  SLOR  method. 

(b)  The  Time-Integration  Approach 


As  a further  simplification  of  Equation  1,  the  frequency  of  the  transonic 

n,  2 ^ 2/3 

flow  can  be  assumed  as  low  so  that  kc  ^ 1 - ^ t ' <<  1 . Equation  1 may 

then  be  reduced  to 

2 k M2  6 . = V <b  + <|>  (3) 

c « yxt  c vxx  Tyy  v ' 

where  t Is  the  thickness-chord  ratio  of  the  airfoil. 

Equation  3 is  suitable  for  the  time-integration  approach.  This  approach 

Is  based  on  the  finite-difference  scheme  that  integrates  Equation  3 in  time 

for  harmonic  aerodynamic  motions  until  the  transient  states  In  the  solution 

disappear  and  the  forces  become  periodic. 


This  method  can  handle  mixed  flows  and  shocks  within  the  accuracy  of 
the  finite-difference  mesh  (Reference  15).  Computations  of  flows  about  two- 
dimensional  airfoils  oscillating  at  supercritical  Mach  numbers  were  performed 
by  Magnus  and  Yoshlhara  (Reference  16)  and  Beam  and  Warming  (Reference  17). 
Such  time  dependent  calculations  Involve  no  Inherent  limiting  assumption  so 
they  are  quite  expensive  In  computing  time. 

A semi-implicit  scheme  was  used  by  Ballhaus  and  Lomax  (Reference  18)  to 
obtain  the  time-accurate  solution  to  the  low-frequency  unsteady  small -distur- 
bance transonic  flow  potential  equation.  Fully  implicit  schemes  were  devel- 
oped by  Beam  and  Warming  (Reference  19)  and  Ballhaus  and  Steger  (Reference  20) 
which  permit  time-step  selection  based  on  accuracy  rather  than  stability 
considerations.  More  recently,  a computer  code,  LTRAN2,  was  developed  by 
Ballhaus  and  Goorjlan  (Reference  21)  which  solves  the  two-dimensional, 
nonlinear,  low-frequency,  small -disturbance  transonic  flow  equation  by  an 
alternating-direction  implicit  algorithm.  The  program  showed  considerable 
Improvement  In  computing  efficiency. 

A time-marching  semi-implicit  finite-difference  technique  was  used  by 
Isogal  (Reference  22)  to  solve  the  full  potential  equation  for  unsteady  tran- 
sonic flow  fields  around  oscillating  airfoils.  The  shock-wave  motions  were 
properly  captured  using  the  quasi -conservative  scheme  proposed  by  Jameson 
(Reference  23)  and  Bauer  and  Korn  (Reference  24)  for  the  steady  flow  problem. 
Numerical  examples  showed  that  for  weak  shock  waves,  the  method  could  give 
a reliable  result  for  a thick,  blunt-nosed  airfoil,  such  as  a modern  super- 
critical airfoil 

In  a recent  study  (Reference  25),  Yu,  Seebass,  and  Ballhaus  improved  the 
LTRAN2  program  by  Including  an  implicit  shock-fitting  algorithm  and  an 
analytically  stretched  coordinate  system.  Shock  waves  were  treated  as 
discontinuities  normal  to  the  free  stream.  Results  were  given  for  an  NACA 


64A006  airfoil  with  oscillating  quarter-chord  flap.  In  Reference  26 
Fung,  Yu,  and  Seebass  provided  a method  that  accounted  for  shock  wave 
motions  due  to  arbitrary  but  unsteady  changes  in  the  boundary  conditions. 

Thus  both  harmonic  and  indiclal  responses  may  be  determined.  Similar 
work  was  also  carried  out  by  Nixon  (Reference  27). 

(c)  Indicia!  Method 

In  solving  for  the  flutter  eigenvalue  equations.  It  is  required  to 
evaluate  the  aerodynamic  coefficients  for  a wide  range  of  reduced  frequencies 
and  different  modes  of  motion  such  as  plunging,  pitching,  and  aileron  pitching, 
etc.  For  each  given  Mach  number,  reduced  frequency,  and  mode  of  motion, 
both  the  aforementioned  relaxation  and  time-integration  methods  require  the 
solution  of  complete  flow  fields  in  order  to  provide  the  aerodynamic  coefficients. 
Both  methods  are  thus  highly  time-consuming  for  the  purpose  of  flutter 
analysis. 

For  the  purpose  of  generating  aerodynamic  coefficients,  a computationally 
more  efficient  approach,  the  indicial  method,  was  used  by  Ballhaus  and 
Goorjian  (Reference  28).  This  method  was  commonly  used  in  the  case  of 
subsonic  flow  (see,  for  example.  Reference  29).  In  the  transonic  application, 
the  method  has  a major  advantage  that  the  aerodynamic  coefficients  can  be 
computed  for  a given  Mach  number,  mode  of  motion  and  a wide  range  of  reduced 
frequencies  by  performing  only  a single  flow  field  computation,  the  indicial 
response. 

The  indicial  response  is  the  flow-field  response  to  a step  change  in  the 
given  mode  of  motion,  and  It  is  computed  using  a time-accurate  finite-difference 
scheme.  From  the  indicial  response,  the  solution  for  all  oscillation  frequencies 
can  be  obtained  by  computing  the  simple  Integrals  in  the  form  of  Duhamel's 
Integral  at  a small  fraction  of  the  cost  of  the  complete  flow-field  solution 
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needed  by  the  relaxation  approach.  For  the  Indlclal  approach,  the  unsteady 
solution  must  be  treated  as  a small  linear  perturbation  about  some  nonlinear 
steady-state  solution. 

The  Indlclal  approach  was  described  In  Reference  28  by  Ballhaus  and 
Goorjlan  and  applied  for  an  aeroelastlc  response  analysis  of  an  NACA  64A006 
airfoil  pitching  In  the  small -disturbance  unsteady  transonic  flow. 

(d)  Other  Methods 

Another  promising  approach  for  both  the  steady  and  unsteady  transonic 
aerodynamics,  which  Is  In  a separate  category.  Is  the  finite  element  approach. 
The  development  of  this  approach  Is  still  In  the  early  stage.  An  account  of 
such  developments  can  be  found  In  a special  section  In  Reference  30  and 
Reference  31. 

In  a recent  report,  Tljdeman  (Reference  32)  presented  wind  tunnel  results 
of  a conventional  airfoil  with  an  oscillating  flap  and  a supercritical  airfoil 
oscillating  In  pitch  in  high  subsonic  and  transonic  flow.  The  current  status 
of  unsteady  transonic  flow  theory  was  reviewed  and  his  test  data  were  used 
to  evaluate  some  of  the  recently  developed  computational  methods. 
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SECTION  III 


AEROELASTIC  EQUATIONS  OF  MOTION  FOR  TWO-DIMENSIONAL 
AIRFOIL  WITH  PITCH  AND  PLUNGE  DEGREES  OF  FREEDOM 

In  the  present  flutter  analysis  of  two-dimensional  airfoils  in  the 
transonic  flow,  it  is  assumed  that  the  airfoils  oscillate  with  two  de- 
grees of  freedom,  plunge  and  pitch.  The  aeroelastic  equations  of  motion 
will  be  derived  in  this  section. 

Figure  1 shows  an  airfoil  plunging  with  a positive  displacement  h (down) 

and  pitching  with  a positive  angle  a (with  the  nose  up)  about  the  elastic 

axis.  The  elastic  axis  is  located  at  a distance  a.b  from  the  midchord  while 

h 

the  mass  center  is  located  at  a distance  x b from  the  elastic  axis.  Both 

a 

distances  are  defined  as  positive  when  measured  toward  the  trailing  edge.  Fi 
example,  if  the  elastic  axis  is  at  1/4  chord  from  the  leading  edge,  ah  = -0.5. 

The  aeroelastic  equations  of  motion  are  derived  by  setting  up  the 
equilibrium  equations  for  three  physical  quantities:  the  inertia  forces, 
the  elastic  forces,  and  the  aerodynamic  forces. 

The  total  vertical  inertia  force  per  unit  span  is  equal  to 

- / (ii  + ra)dm  = - (mh  + Soi)  (4) 

where 

m = total  mass  of  the  airfoil  with  unit  span; 
r = distance  from  elastic  axis,  positive  when  measured  toward 
trailing  edge 

dm  = a differential  element  of  mass;  and 
S = / rdm  = static  moment  of  the  airfoil  about  the  elastic  axis. 
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The  vertical  elastic  force  acting  at  the  airfoil  Is  obtained  as  the 
product  of  the  plunging  displacement  h and  the  bending  stiffness  of  the 
cantilever  wing.  Such  bending  stiffness  Is  customarily  assumed  as  a 
linear  elastic  spring  constant  K^.  Thus  the  elastic  restoring  force  is 
-hK.  . 

n * . 

If  the  aerodynamic  force  is  assumed  as  positive  down  (in  the 
positive  h direction), the  equilibrium  equation  for  all  the  vertical  forces 
is 

mh  + Sa  + hKh  = Qh  (5) 

The  total  inertia  moment  per  unit  span  due  to  pitching  about  the 
elastic  axis  is 

- / r ( h + ra)dm  = - (I^oi  + Sh)  (6) 

p 

where  I = / r dm  = mass  moment  of  inertia  of  the  airfoil  about  the 
a J 

elastic  axis. 

The  elastic  torsional  moment  acting  at  the  airfoil  is  obtained  as 
the  product  of  the  angle  of  pitch  and  the  torsional  stiffness  of  the 
cantilever  wing.  Such  torsional  stiffness  is  assumed  as  a linear  elastic 
spring  constant  Ka>  Thus  the  elastic  restoring  moment  is  - ccK. 

If  the  aerodynamic  moment  Qa  acting  at  the  elastic  axis  is  assumed 
as  positive  in  the  positive  a direction,  the  equilibrium  equation  for  all 
the  moments  is 

Sh  + I a + oK  = Q (7) 

a a a 

Equation  5 can  be  written  in  a nondimensional  form  by  dividing  by  mb. 
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where  E * h/b  Is  the  nondlmenslonal  plunging  displacement;  u^  * (Kh/m)'s  15 
the  uncoupled  plunging  frequency;  and  xq  = S/mb  Is  a coefficient  to  be 
multiplied  by  the  semichord  length  to  locate  the  mass  center  from  elastic 
axis. 

Assuming  harmonic  oscillation,  E and  a may  be  written  as 


5 =C0e 


iut 


a = a e 
0 


lut 


(9) 


where  u Is  the  circular  frequency  of  oscillation  and  E and  a are  the 

oo 

amplitudes  In  plunge  and  pitch,  respectively.  Following  the  definitions 
of  Equation  9,  Equation  8 becomes 

• Co  ‘ Vo  + co(a)h/a,)2  = Qhe’1u>t/mbu2  (10) 

For  more  general  application,  the  right  hand  side  of  Equation  10 
Is  rewritten  to  Include  the  customary  aeroelastlc  constant  p,  q,  and  k^, 

" «o  ‘ Vo  + 5o(a,h/a,)2  = V"ia,t/27rb^kb2  (11) 

2 2 
where  y = m/irb  p Is  the  airfoil  - air  mass  ratio;  q = pU  /2  Is  the 

dynamic  pressure;  and  k^  = ub/U  Is  the  reduced  frequency. 

The  aerodynamic  lifting  force  Qh  is  attributed  to  both  plunging  and 

pitching  motions.  It  can  be  obtained  hy  superposing  the  effects  of  the 

two  motions  through  the  use  of  aerodynamic  coefficients  c£6  and  c^. 

Such  superposition  is  valid  when  the  aerodynamic  equations  considered 

are  linear.  When  the  aerodynamic  equations  considered  are  nonlinear,  such 

superposition  may,  however,  still  be  valid  under  the  assumption  of 


I 


linearization  of  the  unsteady  part  of  the  flow.  The  linearization  Is 
valid  when  the  displacements  are  assumed  to  be  small.  Such  assumption  Is 
In  agreement  with  the  small  perturbation  aerodynamic  theory  used  to  solve 
the  present  nonlinear  transonic  flow  equations. 

The  lifting  force  Qb  can  then  be  written  as 

"h  ■ - <C«S  * cla“>‘'c  (12> 

where  c„.  is  the  lifting  coefficient  of  the  airfoil  due  to  plunging; 

£6 

c is  the  lifting  coefficient  of  the  airfoil  due  to  pitching  about  the 
elastic  axis;  and  6 = h/c  = ?/2.  It  should  be  noted  that  the  lifting 
force  is  defined  as  positive  upward  in  the  aerodynamic  equations  and  this 
sign  convention  is  used  in  the  presentation  of  the  subsequent  tables  of 
aerodynamic  coefficients.  On  the  other  hand,  the  lifting  force  is  defined 
as  positive  downward  in  the  present  flutter  equations.  Thus  there  is  a 
negative  sign  in  Equation  12. 

Assuming  that  the  airfoil  oscillates  harmonically  with  circular  fre- 
quency w yields 

■ - *1"t  < W2  + Vo>«c  (13> 

Substituting  Equation  13  into  Equation  11  and  rearranging  the  terms  gives 

ukb2(co  + Vo}  ‘ 7 (c£6eo/Z  + C£aao) 

= Vkb2(a>h/a))2  (14) 

A similar  procedure  can  be  used  to  treat  the  second  equation  of  mo- 

2 

tion.  Dividing  Equation  7 by  mb  yields 

x c + rZa  + rZwZa  = Q /mb^  (15) 

a a a a a 


I 


t 
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where  r ■ (I  /mb2)3*  Is  the  coefficient  to  be  multiplied  by  the  semi- 


a  a 


chord  to  give  the  radius  of  gyration  about  the  elastic  axis;  and 
id  = ( K /I  V Is  the  uncoupled  pitching  frequency. 


a a a 


Assuming  harmonic  oscillation  gives 


■ Xo5o  ■ *0  + ra(ua/a))2“o  = Qae“iait/mb2u)2. 
= 2 Qae"1a)t/Tiukb2qc2 


Based  on  the  aforementioned  assumption  of  linear  superposition, 
the  aerodynamic  moment  can  be  obtained  as 


2 2 

Q = qc  cm.<5  + qc  Cct 
a mo  ma 


where  c . and  c are  the  moment  coefficients  of  the  airfoil  about  the 
mo  ma 

elastic  axis  due  to  plunging  and  pitching  (about  the  elastic  axis), 
respectively. 

For  harmonic  oscillation 

% ■ o1"*' Vo'2  ♦ V.1'1'  <’«> 

Substituting  Equation  18  Into  Equation  16  and  rearranging  the  terms  gives 

wk.2(x  £_  + r2a  ) + \ [c _oCrt/2  + c a ) 
b oo  oo  ir  mo  o in  o 


r2(to  /“)2o.ykK2 
a a o b 


Grouping  the  terms  In  Equations  14  and  19,  the  final  eigenvalue 
equations  for  flutter  analysis  of  an  airfoil  with  plunging  and  pitching 
degrees  of  freedom  are  obtained. 


Ek>]  - [a]]l0}  ■ MM 
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where  [M],  [A],  and  [k]  are  termed  here  as  mass,  aerodynamic,  and  stiff- 
ness matrices,  respectively,  and  they  are  defined  as 


where  u>  Is  a reference  frequency, 
r 

The  eigenvalue  X Is  a complex  number  defined  as 


1 

X 

a 

2 

X 

r 

a 

a 

where  g Is  the  so-called  structural  damping  coefficient  which  is  assumed 

I 

to  be  small  (as  defined  In  Reference  33  ) and  It  is  assumed  to  be  the  same 

i L 


for  both  plunging  and  pitching  modes. 


SECTION  IV 


r. 


* 

!( 


1 


It 


i'l 


TRANSFORMATION  OF  AERODYNAMIC  COEFFICIENTS 


In  obtaining  the  aerodynamic  coefficients.  It  Is  a common  practice 
to  pitch  the  airfoil  about  the  elastic  axis  so  that  such  coefficients 
can  be  directly  used  to  compute  the  aerodynamic  matrix  [A]  for  aeroelastlc 
analysis.  It  is,  however,  not  uncommon  that  the  pitching  axis  for  ob- 
taining the  aerodynamic  coefficients  is  different  from  the  elastic  axis. 

In  that  case,  the  aerodynamic  coefficients  with  reference  to  the  pitch- 
ing axis  must  be  transformed  to  those  with  reference  to  the  elastic  axis 
for  aeroelastic  analysis.  It  is  the  purpose  of  this  section  to  derive 
a set  of  transformation  equations. 

The  following  derivations  of  the  transformation  equations  for  aero- 
dynamic coefficients  between  two  points  of  pitching  are  valid  only  for 
linear  aerodynamic  theory.  In  the  transonic  region  where  aerodynamic 
equations  are  nonlinear,  the  present  transformation  equations  are  valid 
only  under  the  assumption  of  small  linear  disturbances  about  a nonlinear 
steady  flow.  Also  it  may  be  noted  that  the  effect  of  local  angle  of 
attack  which  depends  on  the  location  of  the  point  of  pitching  of  the  airfoil 
for  a nonlinear  case  is  not  accounted  for  in  the  present  derivation  of 
transformation  equations. 

Figure  2 shows  an  airfoil  pitching  about  point  01  which  is  at  a dis- 


tance bXp  aft  midchord  point  C.  It  is  assumed  that  the  aerodynamic  co- 


efficients have  been  computed  for  this  pitching  axis  with  the  lift  acting  posi- 
tive downward.  The  aerodynamic  lifting  force  and  moment  at  point  0^  are 
designated  as  and  Qa1 , respectively.  The  lifting  and  moment  aerodynamic 
coefficients  at  0^  due  to  plunging  and  pitching  at  0^  are  designated  as  cfc6j, 
c«.al  ’ cm5l*  and  Snal’  resPectively-  ^ desired  to  transform  these  aerodynamic 
coefficients  to  those  at  point  02,  which  is  the  pitching  point  assumed  in 

19 


il 


. -W9  « I , ■ w ■ ^ . 


the  aeroelastlc  analysis  and  It  Is  located  at  a distance  bah  aft  the  mid- 
chord  point  C. 

After  the  transformation,  the  aerodynamic  lifting  force,  moment,  and 

the  four  aerodynamic  coefficients  at  point  02  are  designated  as  , 

Q » . c,  , c . , and  c , respectively. 

if>2  t-a2  m6^  ma2 

Considering  the  equilibrium  of  the  aerodynamic  lifting  forces  and 
moments  at  point  02  gives 


<ah  - V<c'2>\  ♦ y 


(23) 


I 

I 


Defining  that 
s = (ah  “ xp)/2 

gives 


For  small  angle  of  pitch  and  rigid  airfoil 

h2  = h.|  + sca^ 
a2  = a, 


(24) 


(25) 


(26) 


With  the  definitions  that 


Equations  26  can  be  written  in  a nondimenslonal  matrix  form 


^ ^ ■ : 


Substituting  Equations  30  and  32  Into  the  equilibrium  equation  25 
yields 


1 

1 0 

sc 

| 1 

00 

CM 

for  i 

6^  and  , and  then  rearranging  the  equation  gives 

rp]  = _J — [.  °. , fp-|  L s 
LN  sc  1 LtJ  0 1 


r pi  = 

1 : 

0 

i 

L*  J 

sc  | 

i 

1 

rfn 

fl  2 

I 

! 

CM 

4- 

1 

f22_ 

where 


fn  ■ 

f12  * - "“ta, 

f21  ’ * "c2cmSl 

f22  = ’ “I c2scta,  ' 

♦ <ic2c 


Thus  matrix  [F]  is  defined  In  both  Equations  33  and  36.  Equating  the 
corresponding  entries  In  the  two  matrices  gives 


Cl&2 

* c 

t61 

c. 

* c 

- sc 

vla2 

t61 

cm«2 

s r 

m«1 

- sclfil 

c 

* c 

- s(c. 

ma2 

ma  i 

'la 

The  terms  sc„r  , sc  . , and  s ca.  are  relatively  small  as  compared 


A/\j  ^ mu  i 

with  other  terms.  If  these  three  terms  are  neglected.  Equations  37  are 


reduced  to  precisely  the  same  form  as  that  presented  in  Reference  5 by 


Traci  et  al 


SECTION  V 

FLUTTER  SOLUTION  PROCEDURE 

Methods  for  the  solution  of  flutter  eigenvalue  equations  are  avail- 
able in  common  texts  of  aeroelasticity  (see,  for  example.  Reference  33). 

In  this  study,  the  standard  U-g  method  is  used.  The  basic  idea  of  the 
U-g  method  is  to  search  for  a set  of  aeroelastic  parameters  that  result 
in  a positive  real  eigenvalue  x or  a zero  structural  damping  coefficient 
9- 

In  the  solution  of  the  present  flutter  eigenvalue  Equation  20,  the 
mass  matrix  [M]  and  the  stiffness  matrix  [K]  can  be  obtained  for  a set  of 
specified  parameters:  airfoil-air  mass  ratio  y;  location  of  the  mass  cen- 
ter x^;  radius  of  gyration  ra;  plunging  frequency  ay;  pitching  frequency 
ay;  and  a reference  frequency  ay  (used  for  nondimensionalization).  In  this 
study,  ay  is  always  assumed  as  unity. 

To  obtain  the  aerodynamic  matrix  [A],  aerodynamic  coeffi cents  c^, 

cia’  cm<5’  an(1  cmo  must  be  first  comPuted  based  on  the  following  given 

conditions:  airfoil  configuration;  free  stream  Mach  number  M;  location 

of  the  point  of  pitching  x , and  reduced  frequency  k . If  the  point  of 

P c 

pitching  does  not  coincide  with  the  elastic  axis,  the  four  aerodynamic 
coefficients  must  be  transformed  to  the  elastic  axis.  Such  transformation 
procedure  has  been  given  in  detail  in  Section  IV. 

After  the  three  matrices  [M],  [A],  and  [K]  are  found.  Equation  20 
can  be  solved  directly  and  the  two  complex  eigenvalues  and  x 2 can  be 
obtained.  The  flutter  phenomenon  is  obtained  when  either  of  the  two 
eigenvalues  is  positive  and  real  or  when  either  of  the  two  corresponding 
structural  damping  coefficients  g is  zero.  The  eigenvector  corresponding 


a 


ii 


I 


i 

\ 


i 
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to  A1  or  Ag,  respectively,  gives  the  ratio  between  plunging  and  pitching 

mode  amplitudes  which  describes  the  flutter  mode. 

In  order  to  search  for  the  flutter  solution,  the  reduced  frequency  kc 

is  varied  in  small  increments,  the  corresponding  aerodynamic  coefficients 

are  computed,  and  the  corresponding  eigenvalues  A's  are  obtained  from 

Equation  20.  The  flutter  solution  is  obtained  when  a positive  real  eigenvalue 

A is  found.  The  repetitive  computations  of  such  transonic  aerodynamic 

coefficients  are  enormously  laborious  and  expensive.  It  appears  desirable 

to  compute  a table  of  aerodynamic  coefficients  based  on  selected  values 

of  k with  appropriate  intervals  so  that  the  aerodynamic  coefficients  for 
c 

any  k value  within  the  acceptable  range  can  be  obtained  by  the  method  of 
interpolation  or  extrapolation.  In  this  study,  typical  values  chosen 
for  kc  are  0.0,  0.05,  0.10,  0.15,  0.2,  0.25,  and  0.3.  The  method  of 
interpolation  is  based  on  the  Lagrange  interpolation  formula  as  described 
in  coirmon  texts  on  numerical  methods  (see,  for  example.  Reference  37). 

In  this  study,  the  eigenvalue  is  defined  as 

A = y(l  + ig)u,2b2/U2  (38) 

Thus  the  search  for  a positive  real  eigenvalue  A becomes  the  search  for 
zero  value  of  g,  the  so-called  structural  damping  coefficient.  As  the  reduced 
frequency  k£  varies  in  small  increment,  the  corresponding  g value  changes. 

The  zero  g value  is  obtained  when  it  changes  sign.  Physically,  the  airfoil 
goes  from  a stable  state  to  an  unstable  state  when  the  g value  changes  from 
negative  to  positive,  and  vice  versa. 

At  the  flutter  points,  g = 0 and 

A * uw2b2/U2  (39) 
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For  general  application,  the  results  for  the  flutter 
frequency  are  nondimenslonalized  as 

U*  * U/bu) 

a 

and 

U)*  = 0 j/lD 

a 

respectively. 

Incorporating  Equation  39  into  Equations  40  and  41  i 
manner  and  using  the  definition  that  kb  = iob/U  gives 


speed  and  flutter 

(40) 


(41) 


an  appropriate 


SECTION  VI 


FLUTTER  RESULTS  FOR  NACA  64A006  AIRFOIL 

(a)  NACA  64A006  Airfoil  Configuration 

In  order  to  perform  accurate  aerodynamic  computations,  it  is  necessary 
to  use  accurate  data  for  the  airfoil  configurations,  especially  in  the  region 
of  the  leading  edge.  The  equations  for  accurate  curve  fitting  and  tabulated 
data  for  eight  AGARD/SMP  standard  airfoil  configurations  were  provided  by 
Dr.  J.J.  Olsen  of  the  Air  Force  Flight  Dynamics  Laboratory  (Reference  38). 

The  equations  for  the  eight  airfoil  configurations  were  given  by  Dr. 
Olsen  as 

Zu(x)  = c ( x ) + 0.5  t ( x ) 

z,<„  - «(*)  - 0.5  t(xi  <«, 

t ( x ) = aQ/x  + a]  + a2x  + a3x  + a4x  + a5x 
c ( x ) = b]  + b2x  + b3x2  + b4x3  + bgx4 

. 

where  Z (x)  and  Z (x)  are  the  equations  for  the  upper  and  lower  surfaces, 

U x» 

respectively;  x is  the  nondimensional  axis  with  leading  edge  at  x = 0 and 
trailing  edge  at  x = 1;  c(x)  is  the  function  of  camber;  t(x)  is  the  function 
of  thickness;  and  constants  a's  and  b's  for  the  eight  airfoils  are  available 
in  Reference  38. 

In  this  study,  an  NACA  64A006  airfoil  was  first  selected  for  flutter 
analysis.  There  is  no  camber  for  this  airfoil.  The  constants  for  the  thick- 
ness function  are:  aQ  = 0.1369014;  a.|  = 0;  a2  = -0.0496509;  a3  = 0.0609015; 
a4  = -0.3455406;  a^  = 0.1976486.  The  airfoil  configuration  was  plotted  in 
Figure  3.  Part  IV  of  Reference  38  contains  revised  coefficients,  however  the 

above  values  are  suitable  for  purpose  of  illustration. 
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NfiCO  64P006 


Figure  3.  NACA  64A006  Airfoil  Configuration 


1 

(b)  Aerodynamic  Coefficients 

In  this  study,  the  airfoil  was  assumed  to  oscillate  with  only  two  degrees 
of  freedom,  plunge  and  pitch.  Consequently,  only  four  aerodynamic  coefficients 
were  computed:  the  lifting  force  due  to  unit  plunging  displacement  c£(S;  the 
lifting  force  due  to  unit  pitching  rotation  c^  ; the  moment  due  to  unit  plunging 
displacement  cm5;  and  the  moment  due  to  unit  pitching  rotation  c . 

The  aerodynamic  coefficients  were  computed  by  using  two  separate  computer 
programs: 

(1)  STRANS2  and  UTRANS2  based  on  the  harmonic  or  relaxation  method 
(Reference  3-8); 

(2)  LTRAN2  based  on  the  indicial  and  time-integration  methods  (Reference  28). 

The  airfoil  was  assumed  to  pitch  about  both  the  quarter-chord  axis  and 

■ 

the  leading  edge.  Free  stream  Mach  numbers  considered  were  0.7,  0.8,  0.85, 

and  0.9,  respectively.  LTRAN2  was,  however,  not  applicable  when  M = 0.9. 

00 

Reduced  frequencies  considered  were,  in  general,  between  0.0  and  0.25  with  an 
increment  of  0.05. 

It  was  suggested  by  Dr.  S.R.  Bland  of  NASA-Langley  Research  Center  that 
for  a lower  free  stream  Mach  number  such  as  0.7,  the  aerodynamic  solutions 
computed  by  using  the  transonic  computer  programs  STRANS2,  UTRANS2,  and  LTRAN2 


* ! I 

i J 


I 

i 


If: 


1 /, 

; 


i, 


should  agree  with  the  aerodynamic  solution  computed  by  using  the  linear  flat 
plate  theory  for  compressible  flow.  A copy  of  the  computer  program  based  on 
the  Kernel  function  method,  to  solve  for  the  linear  flat  plate  aerodynamic 
equation,  was  provided  by  Dr.  Sland  of  NASA-Langley. 

Attempts  were  made  to  use  program  LTRAN2  to  analyze  flat  plates  but  were 
unsuccessful.  As  an  alternative,  a thin  parabolic  arc  airfoil  with  t * 0.01 
was  used  to  simulate  a flat  plate.  LTRAN2  has  an  option  (by  setting  y = -1 ) 
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for  the  solution  of  only  the  linear  aerodynamic  equation.  This  option  was 
used.  Although  STRANS2  and  UTRANS2  can  accurately  treat  a flat  plate,  the 
same  thin  parabolic  arc  airfoil  was  chosen  for  the  analysis  so  that  the  results 
obtained  from  both  the  time-integration  program  and  the  relaxation  program 
can  be  compared  directly. 

Figure  4 shows  the  distributions  of  steady  pressure  coefficients  for 
the  thin  parabolic  arc  airfoil  with  x = 0.01  at  M * 0.7  obtained  by  using  both 
STRANS2  and  LTRAN2.  The  two  curves  are  in  good  agreement.  It  is  noted  that 
for  the  steady  solution,  both  programs  STRANS2  and  LTRAN2  are  based  on  the  same 
relaxation  method. 

Table  1 shows  the  aerodynamic  coefficients  for  the  thin  parabolic  arc 
airfoil  with  x = 0.01  pitching  about  the  quarter  chord  axis  at  M = 0.7  obtained 
by  using  both  the  relaxation  method  (UTRANS2)  and  the  indicial  method  (LTRAN2). 
A linear  flat  plate  solution  obtained  by  using  the  Kernel  function  method  is 
also  shown  for  comparison.  In  general,  the  results  of  the  relaxation  method 
agree  with  the  linear  flat  plate  solution  better  than  the  results  of  the  indi- 
cial method.  This  can  be  attributed  to  two  reasons.  First,  a relatively 
coarse  mesh  (35  x 38)  was  used  in  the  indicial  method.  Second,  the  low  fre- 
quency approximation  (by  omitting  the  <|>tt  term  in  the  aerodynamic  equation) 
used  in  LTRAN2  is  less  accurate  when  the  free  stream  Mach  number  deviates 
farther  from  Mw=  1.0.  The  term  is,  however,  retained  in  the  relaxation 
program  UTRANS2. 

Figures  5 and  6 show,  respectively,  the  real  and  imaginary  parts  of  the 
distributions  of  the  difference  in  unsteady  pressure  coefficients  across  the 
airfoil.  It  Is  seen  that  all  the  three  solutions  obtained  from  relaxation 
method,  time-integration  method,  and  Kernel  function  method,  respectively. 
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Figure  4.  Distribution  of  Steady  Pressure  Coefficients  for 
Thin  Parabolic  Airfoil  (t  = 0.01)  at  M = 0.7 


TABLE  1 

Aerodynamic  Coefficients  for  flat  plates  pitching  at  the  quarter  chord  at 
M = 0.7 


Reduced  Frequency 


.0 

kc=< 

1.05 

kc=o.io 

k =0.15 
c 

k =0.20 
c 

k *0.25 
c 

Imag 

Real 

Imag. 

Real 

Imag. 

Real 

Imag. 

Real 

Imag. 

Real 

Imag. 

0.0 

0.059 

0.462 

0.188 

0.868 

0.330 

1.220 

0.532 

1 .540 

0.693 

1 .780 

0.0 

0.055 

0.385 

0.153 

0.712 

0.271 

0.999 

0.425 

1.210 

0.422 

1 .510 

0.0 

0.062 

0.397 

0.167 

0.719 

0.272 

0.988 

0.365 

1 .222 

0.442 

1 .434 

0.0 

9.230 

-1.170 

8.680 

-1.880 

8.140 

-2.200 

7.680 

-2.660 

7.120 

-2.770 

0.0 

8.078 

-0.382 

7.212 

-1.204 

6.715 

-1 .445 

6.294 

-1 .603 

5.778 

-1 .634 

y 

7.980 

-1.033 

7.282 

-1.308 

6.739 

-1 .323 

6.323 

-1.222 

6.004 

-1.066 

EH 

-0.002 

-0.021 

-0.004 

-0.039 

-0.006 

-0.057 

-0.010 

-0.075 

-0.012 

-0.090 

0.0 

0.002 

-0.009 

0.006 

-0.021 

0.013 

-0.032 

0.021 

-0.045 

0.042 

-0.056 

0.0 

0.002 

0.000 

0.008 

-0.002 

0.019 

-0.004 

0.032 

-0.008 

0.048 

-0.013 

0.0 

-0.409 

.029 

-0.394 

0.042 

-0.379 

0.041 

-0.373 

0.049 

-0.358 

0.047 

0.0 

-0.207 

-0.070 

-0.204 

-0.114 

-0.211 

-0.168 

-0.223 

-0.219 

-0.231 

-0.261 

0.0 

-0.005 

-0.067 

-0.014 

-0.140 

-0.023 

-0.208 

-0.031 

-0.273 

-0.038 

-0.337 

Method  1:  Parabolic  arc,  t = 0.01,  Indicial  method  (35  x 38) 
Method  2:  Parabolic  arc,  t * 0.01,  Relaxation  method  (59  x 60) 
Method  3:  Linear  flat  plate  theory.  Kernel  function  method. 
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Figure  5.  Distribution  of  Unsteady  Pressure  Coefficients  (Real 
Part)  for  Plates  Pitching  about  the  Quarter  Chord 
with  kh  - 0.05  and  M - 0.7. 
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Figure  6.  Distribution  of  Unsteady  Pressure  Coefficients 
(Imaginary  Part)  for  Plates  Pitching  about  the 
Quarter  Chord  with  k.  = 0.05  and  M * 0.7 


are  In  good  agreement. 

Two  sets  of  aerodynamic  coefficients  for  an  NACA  64A006  airfoil  pitching 
about  the  quarter  chord  axis  at  M = 0.7,  0.8  and  0.85  are  given  In  Tables  2,  3, 
and  4,  respectively,  for  various  values  of  reduced  frequency  kc.  One  set  was 
obtained  by  using  the  relaxation  programs  STRANS2  and  UTRANS2  and  the  other 
set  by  the  indicial  program  LTRAN2.  For  M = 0.9,  only  results  by  the  relaxa- 
tion programs  were  obtainable  and  they  are  given  in  Table  5.  It  is  difficult 
to  draw  conclusions  consistent  among  all  tables  on  the  comparison  of  the  two 
sets  of  results.  However,  the  agreements  in  real  parts  are,  in  general,  better 
than  those  in  imaginary  parts.  The  agreements  in  lifting  coefficients  are, 
in  general,  better  than  those  in  moment  coefficients. 

Figures  7,  8,  9,  and  10  show  the  distributions  of  steady  pressure  coeffi- 
cients across  the  airfoil  for  the  free  stream  Mach  numbers  of  0.7,  0.8,  0.85, 
and  0.9,  respectively.  In  both  programs  STRANS2  and  LTRAN2,  the  steady  solu- 
tions are  based  on  the  relaxation  method.  The  small  difference  between  the 
two  steady  pressure  curves  in  Figs.  7 and  8 may  be  due  to  the  difference  in 
mesh  size. 

For  the  case  M = 0.7  and  0.8,  the  shock  is  not  present.  It  is  present  at 
M = 0.85  as  shown  in  Figure  9.  A finer  mesh  (79  x 99)  was  thus  used  in  LTRAN2 
for  this  case.  An  experimental  solution  by  Tijdeman  (Reference  32)  was  also 
plotted  for  comparison.  A pronounced  shock  is  observed  in  the  case  M = 0.9 
as  shown  in  Figure  10. 

It  is  noted  that  a shock-point-operator,  included  in  STRANS2  to  improve 
accuracy  In  the  calculation  of  strong  shock  waves  (supersonic  to  subsonic 
transition)  for  supersonic  free  stream  flow  was  taken  out  in  the  present 
computation  of  the  steady  pressure  distributions.  As  pointed  out  by  Holman 
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TABLE  2 

Aerodynamic  Coefficients  for  NACA  64A006  pitching  at  the  quarter 
chord  at  M = 0.7 


Reduced  Frequency 


5 

kc  = 0.15 

HI 

.20 

HO 

Imaq.  Real  Ima 


0.0  0.061  0.474  0.196  0.884  0.351  1.260  0.552  1.540  0.693  1.820 

0.0  ' 0.0  0.043  0.438  0.165  0.759  0.288  1.056  0.468  1.283  0.628  1.466 


10.000  0.0  9.470  -1.220  8.840  -1.960  8.390  -2.340  7.710  -2.760  7.270  -2.7 

8.845  0.0  8.693  -0.547  7.769  -1.403  7.073  -1.665  6.545  -1.845  5.983  -1.8 


0.0  0.0  -0.001  -0.021  -0.004  -0.040  -0.005  -0.059  -0.008  -0.076  -0.008 
0.0  0.0  0.002  -0.011  0.006  -0.021  0.014  -0.033  0.028  -0.050  0.032 


-0.435  0.0  -0.417  0.026  -0.401  0.036  -0.394  0.033  -0.382  0.040  -0.372  0.033 
-0.217  0.0  -0.217  -0.071  -0.217  -0.120  -0.219  -0.175  -0.231  -0.226  -0.241  -0.2 


m 


Method  1:  Indicial  Method  (35  x 38) 
Method  2:  Relaxation  Method  (59  x 60) 


TABLE  4 

Aerodynamic  Coefficients  for  NACA  64A006  pitching  at  the  quarter  chord 


k *0.0  k -0.05 
c c 


Reduced  Frequency 


k *0.10  k =0.15 

c c 


k *0.18 
c 


Real 

Imag. 

Real 

Imag. 

Real 

Imag. 

Real 

Imag. 

0.136 

0.666 

0.492 

1.180 

0.938 

1.470 

1.187 

1.525 

0.163 

0.626 

0.480 

0.962 

0.783 

1 .148 

0.832 

1 .039 

13.310 

-2.720 

11.820 

-4.920 

9.770 

-6.250 

8.474 

-6.596 

11.803 

-3.701 

9.160 

-4.491 

7.223 

-4.076 

5.804 

-3.689 

-0.001 

-0.034 

-0.005 

-0.068 

-0.012 

-0.102 

-0.018 

-0.122 

0.004 

-0.026 

0.009 

-0.054 

0.013 

-0.089 

0.002 

-0.096 

-0.676 

0.021 

-0.677 

0.053 

-0.677 

0.082 

-0.677 

0.101 

-0.501 

-0.071 

-0.510 

-0.147 

-0.514 

-0.210 

-0.498 

-0.180 

kc.0.2S 

Real 

Imag. 

1.570 

1.480 

Method  1:  Indlclal  Method  (79  x 99) 
Method  2:  Relaxation  Method  (59  x 60) 





TABLE  5 


Aerodynamic  Coefficients  for  NACA  64A006  pitching  at  the  quarter  chord 
at  M * 0.90 


Reduced  Frequency 


kc=( 

).0 

II 

) . 05 

kc=( 

) . 07 

| k =0.10 

kc=( 

) . 1 2 

Real 

I mag 

Real 

I mag. 

Real 

Imag. 

Real 

Imag. 

Real 

Imag. 

0.0 

0.0 

0.095 

0.371 

0.153 

0.481 

0.247 

0.617 

0.309 

0.681 

8.852 

0.0 

7.444 

-1.757 

6.898 

-1.990 

6.098 

-2.117 

5.527 

-2.047 

0.0 

0.0 

-0.016 

-0.061 

-0.027 

-0.080 

-0.047 

-0.104 

-0.062 

-0.116 

1.437 

0.0 

-1.228 

0.260 

-1.150 

0.299 

-1 .031 

0.334 

-0.940 

0.349 

Relaxation  Method  (59  x 60) 
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Figure  7.  Distribution  of  Steady  Pressure  Coefficients  for 
NACA  64A006  Airfoil  at  M = 0.7 
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Figure  9.  Distribution  of  Steady  Pressure  Coefficients  for 
NACA  64A006  Airfoil  at  M = 0.85 
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Figure  10.  . Distribution  of  Steady  Pressure  Coefficients  for 
NACA  64A006  Airfoil  at  M = 0.9 
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1 


(Reference  8),  the  shock- point-operator  Is  essential  only  for  low  supersonic 
flow.  It  leads  to  a strong  over  expansion  of  the  shock  waves  for  high  sub- 
sonic flows. 

In  Reference  5,  Traci  et  aK  provided  the  aerodynamic  coefficients  for 
an  NACA  64A006  airfoil  pitching  about  the  leading  edge.  Their  results  and 
the  present  results  for  = 0.8,  0.85,  0.9,  and  0.95  are  presented  In  Tables 
6,  7,  8,  and  9,  respectively,  for  comparison. 

Table  6 shows  the  results  at  = 0.8  obtained  by  using  the  indlclal  method, 
the  time- Integration  method  with  two  meshes,  the  relaxation  method,  and  those 


provided  in  Reference  5.  The  results  obtained  by  using  the  time-integration 
method  with  the  finer  mesh  (79  x 99)  should  be  the  most  accurate  ones. 

In  Table  7,  the  results  obtained  for  = 0.85  by  using  the  Indlclal, 
time-integration,  and  the  relaxation  methods  compare  fairly  well  with  those 
provided  in  Reference  5.  In  Tables  8 and  9,  the  results  obtained  for  = 0.9 
and  0.95,  respectively,  by  using  STRANS2  and  UTRANS2  agree  well  with  those 
obtained  in  Reference  5. 

It  is  noted  that  the  aerodynamic  coefficients  obtained  by  using  the  relaxa- 
tion method  drop  off  in  high  ^-values.  This  may  be  due  to  the  fact  that  the 
convergence  becomes  more  difficult  in  higher  k£  values. 

(c)  Flutter  Results 

In  this  study,  the  flutter  speed  is  predicted  by  using  the  standard  U-g 
method  as  described  in  Section  V.  A typical  U-g  plot  is  given  in  Figure  11 
for  the  case  of  an  NACA  64A006  airfoil  pitching  about  the  1/4-chord  axis. 

In  this  flutter  analysis,  the  effects  of  four  parameters  were  considered: 


airfoil-air  mass  ratio  p;  position  of  the  mass  center  xa;  position  of  the 
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TABLE  6 


Aerodynamic  coefficients  for  NACA  64A006  pitching 
at  the  leading  edge  at  M = 0.8 


' Bp 

Reduced  Frequency  Ur. 

) 

Method 

HHi) 

.0 

0.05 

0.10 

0.15 

0.20 

rife 

unnn 

mi 

Imag. 

Real 

1EE1I 

Imag. 

■roil 

Imag. 1 

1 

0.0 

[<M*I 

mmm 

UMLUM 

wmtkA 

mw-M 

Bt  •7*1*1 

i&KUi 

r 

2 

.091 

.574 

.328 

IBSa 

11 

■R51 

\ 

3 

.309 

.951 

6 

4 

.147 

.551 

.259 

.837 

.674 

5 

0.0 

0.0 

.175 

.47 

.259 

.817 

.471 

1.15 

.708 

1 

12.6 

0.0 

Tswm 

-2.52 

9.98 

-3.47 

8.60 

-3.82 

8.00 

-4.13 

r 

11.5 

-1.82 

10.3 

-2.75 

9.04 

-2.94 

8.74 

-2.84 

c* 

9.66 

-2.59 

a 

10.5 

0.0 

9.5 

-1.82 

8.35 

-2.89 

6.58 

-3.42 

. 

H 

11.5 

0.0 

10.0 

-2.08 

8.84 

-2.99 

7.76 

-3.43 

6.94 

-3.55 

1 

0.0 

0.0 

-.027 

-.141 

-.072 

-.256 

-.117 

mm 

-.165 

- .437 

-.047 

-.137 

-.069 

-.257 

-.111 

Bee 

-.167 

- .435 

V 

m6 

msm* 

Pm 

-.246 

fl 

0.0 

0.0 

-.036 

- .159 

-.057 

- .245 

-.146 

- .401 

n 

0.0 

0.0 

-.043 

- .130 

- .055 

- .229 

-.096 

-.332 

-.147 

- .401 

i 

-3.13 

0.0 

-2.82 

.540 

-2.55 

.72 

-2.26 

.78 

-2.18 

.83 

-2.86 

.453 

-2.62 

.37 

-2.37 

.38 

-2.32 

.37 

-2.51 

.40 

a 

■ 

-3.00 

0.0 

-2.72 

.391 

-2.43 

.72 

-2.00 

.83 

mu 

-3.12 

0.0 

-2.75 



.500 

-2.46 

.70 

-2.21 

.79 

-2.02 

.81 

i 

I 


: 


l 


Method  1 - LTRAN2  Indicial  Method  (35  x 38  Grid) 


Method  2 - LTRAN2  Time-Integration  (35x38  Grid) 


Method  3 - 
Method  4 - 
Method  5 - 


LTRAN2  Time- 
STRANS2  and 
STRANS2  and 


Integration  (79  x 99  Grid) 
UTRANS2  (54  x 55  Grid) 
UTRANS2  (Reference  5) 
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TABLE  7 


Aerodynamic  coefficients  for  NACA  64A006  pitching 
at  the  leading  edge  at  M = 0.85 


Real  Ima 


1 0.0 


0. 

05 

Real 

Imag. 

.238 

.689 

.260 

.678 

.184 

.629 

.246 

.575 

■ 


Reduced  Frequency  (k„) 


0.10 


Real  Ima 


C 

>.15 

Real 

Ima 

h. 

.09  .899 

.07  .910 

.969  .715 

.04  .816 


-4.75  10.9 

-4.50  11.4 

-3.50  9.10 

-4.36  10.5 


IliHil 


-3.97 

-3.97 

-3.44 

-3.62 


.13  -3.32 
.06  -3.41 
.823  -2.86 
.04  -3.03 


85 

23  1-2.63 


Method  1 
Method  2 
Method  3 
Method  4 


LTRAN2  Indlclal  Method  (35  x 38  Grid) 
LTRAN2  Time  Integration  (35  x 38  Grid) 
STRANS2  and  UTRANS2  (54  x 55  Grid) 
STRANS2  and  UTRANS2  (Reference  5) 


-5 

.99 

-4 

.56 

-3 

.81 

-5 

.38 

_ 

.455 

- 

.457 

- 

.355 

- 

.356 

1 

.38 

.867 

.73 

1 

.23 

Real 

1 

.22 

1 

.20 

.934 

.941 

7 

.26 

8 

.44 

6 

.29 

7 

.52 

0.20 


Ima 


-.516 

-.522 

-.394 

-.428 


Aerodynamic  Coefficients  for  NACA  64A006  pitching 

{ 

at  the  leading  edge  at  M a 0.95 


Reduced  Frequency  (kc) 

mm 

m 

0 

05 

0.10 

0.15 

Method 

Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

c* 

1 

0.0 

0.032 

0.284 

0.043 

0.584 

<5 

2 

Wxm 

K 

0.057 

0.315 

C* 

1 

6.652 

0.0 

6.076 

-0.180 

a 

2 

7.095 

Bra 

KlETiTij 

cmjt 

1 

m 

I 

-0.018 

-0.255 

5 

2 

0.0 

0.0 

■ H 1 

Cm 

1 

-2.932 

0.0 

0.170 

-2.675 

-0.004 

a 

2 

-3.050 

0.0 

9 

0.640 

Method  1:  STRANS2  and  UTRANS2  (54  x 55  Grid) 
Method  2:  STRANS2  and  UTRANS2  (Reference  5) 


; 


DAMPING  COEFFICIENT 


kcsO.I3 


elastic  axis  «h;  and  the  plunge-pitch  frequency  ratio  (<oh/u>o).  In  all 
the  following  cases,  the  values  for  the  radius  of  gyration  rQ  and  the 
reference  frequency  u were  assumed  as  0.5  and  1.0,  respectively. 


1 . Flutter  Speeds  for  Flat  Plate  Pitching  about 


1 /4-Chord  Axis  at  M 


The  aerodynamic  coefficients  obtained  In  Table  1 were  used  to  perform 
a flutter  analysis  of  the  special  case  of  flat  plates  pitching  about  the 
1/4-chord  axis  at  = 0.7.  The  results  for  flutter  speed  and  the  corresponding 
reduced  frequency  versus  the  airfoil-air  mass  ratio  tf  were  obtained  and 
shown  in  Figure  12.  The  mass  center  was  at  1/8-chord  aft  the  elastic  (pitch) 
axis  and  the  plunge-pitch  frequency  ratio  (u^/a^)  was  0.2. 

Figure  12  shows  that  the  curves  of  flutter  speed  resulted  from  both  the 
indlcial  and  relaxation  methods  are  in  good  agreement.  The  two  curves  are 
higher  than  that  resulted  from  the  Kernel  function  method  by  less  than  10%. 

The  agreement  among  the  three  curves  is  better  in  the  lower  p range.  The 
predicted  flutter  speeds  are,  however,  less  reliable  in  the  lower  p range 
because  the  reduced  frequencies  kc  become  higher  than  0.2. 

It  is  noted  that  the  LTRAN2  results  become  less  accurate  when  k > 0.2. 


Since  the  low  frequency  approximation  was  not  used  in  UTRANS2,  the  valid  k 
range  can  be  extended  somewhat  higher  than  0.2. 


2.  Flutter  Speeds  for  NACA  64A006  Airfoil  Pitching  about  1/4 


Chord  Axis  with  Varying  p at  M = 0.7,  0.8,  and  0.85 


In  this  analysis,  the  effect  of  the  airfoil -air  mass  ratio  u on  flutter 


speed  was  first  studied.  The  range  of  the  p value  was  estimated  by  considering 


two  extreme  cases 


For  an  aluminum  parabolic  arc  airfoil  (t  = 0.06)  with  skin  thickness  of 


0.1%  of  the  chord  length  and  without  stringers  at  sea  level,  the  p value  was 


REDUCED  FREQ,  (k-)  FLUTTER  SPEED  (U/bw) 


Figure  12.  Comparison  of  Flutter  Results  between  Linear  Flat  Plate  Theory 
and  Present  Parabolic  Arc  Approximations  (**  = 0.01)  by  LTRAN2 
and  UTRANS2  at  M * 0.7 


1 
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estimated  as  6.  For  the  same  airfoil  filled  with  solid  aluminum  at  50,000 
feet  above  sea  level,  the  y value  was  estimated  as  760. 

Figure  13  shows  the  curves  for  the  flutter  speed  and  the  corresponding 
reduced  frequency  versus  the  airfoil-air  mass  ratio  y for  an  NACA  64A006 
airfoil  pitching  at  the  1/4-chord  axis  at  = 0.7.  The  curves  were  obtained 
for  two  different  values  of  the  plunge-pitch  frequency  ratios 
0.1  and  0.3.  The  mass  center  was  assumed  to  be  at  1/8-chord  length  aft 
the  elastic  (pitch)  axis.  The  y values  considered  ranged  from  25  to  300. 

Figure  13  shows  that  the  flutter  speed  increases  sharply  with  increasing 
y for  (u>h/u)a)  = 0.1.  Such  increasing  trend  slows  down  rapidly  as  (o^/u^) 
increases.  The  curves  for  flutter  speed  are  lower  for  higher  values  of 

(“h/u,a)- 

It  Is  important  to  mention  that  the  reduced  frequency  kc  increases  as 
y decreases  and  kc  increases  as  (u^/uj^)  increases.  Thus  it  becomes  difficult 
to  obtain  flutter  results  in  the  lower  y region  and  higher  ) values 

since  k goes  out  of  the  valid  range  of  the  present  methods. 

v 

Two  flutter  speed  curves  obtained  by  using  the  linear  flat  plate  theory 
for  (oj^/oj^)  = 0.1  and  0.3,  respectively,  are  also  shown  in  Figure  13.  Each 
curve  is  seen  to  be  lower  but  in  good  agreement  with  the  two  corresponding 
curves  for  the  NACA  64A006  airfoil. 

Using  the  same  conditions  assumed  In  obtaining  Figure  13,  similar 
flutter  results  were  obtained  for  M = 0.8  and  0.85  and  are  shown  in  Figures 

oo 

14  and  15,  respectively.  The  same  general  trends  are  observed.  The  curves 
for  flutter  speed  drop  as  the  free  stream  Mach  number  increases.  The  relaxa- 
tion results  and  the  indicial  results  are  in  good  agreement  in  both  figures. 
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Figure  13.  Effect  of  Mass  Ratio  on  Flutter  Speed  for  NACA  64A006  Airfoil 


for  Various  Frequency  Ratios  at  M = 0.7 


■ 


54 


REDUCED  FREO.(k)  FLUTTER  SPEED(U/b« 


NACA  64A006  , M *0.80 
a.  * -0.5  . x„»  0.25 
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Figure  15.  Effect  of  Mass  Ratio  on  Flutter  Speed  for  NACA  64A006  Airfoil 
for  Various  Frequency  Ratios  at  M = 0.85 
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3,*  Flutter  Speeds  for  NACA  64A006  Airfoil  Pitching  about  1 74- 
Chord  Axis  with  Varying  at  Mm  = 0.7,  0.8,  and  0.85 

* 

In  this  analysis,  the  effect  of  the  position  of  the  mass  center  x on 

a 

flutter  speed  was  studied.  The  NACA  64A006  airfoil  was  assumed  to  pitch  about 
the  1 /4-chord  axis  and  the  airfoil -air  mass  ratio  p was  assumed  as  100. 

Curves  for  flutter  speed  and  the  corresponding  reduced  frequency  versus 
the  location  of  the  mass  center  x^  were  obtained  for  three  different  values 
of  (a^/u^).  The  mass  center  location  was  moved  from  the  1/4-chord  to  1/2- 
chord.  These  curves  are  shown  in  Figures  16,  17,  and  18  for  = 0.7,  0.8, 
and  0.85,  respectively. 

All  three  figures  illustrate  the  same  general  trend.  The  flutter  speed 

decreases  as  the  mass  center  is  moved  from  the  elastic  (pitclt)  axis 

(1/4-chord)  toward  the  1/2-chord.  The  flutter  speed  decreases  as  (w./w  ) in- 

n a 

creases.  The  curves  for  the  reduced  frequency  demonstrate  a behavior  opposite 
to  that  of  the  flutter  curves.  Again,  the  curves  for  flutter  speed  drop 
slightly  as  the  free  stream  Mach  number  increases. 

It  is  important  to  point  out  that  the  agreements  in  flutter  speed  between 
the  indicial  and  the  relaxation  methodsare  good  in  all  three  figures.  However, 
such  agreements  are  better  when  the  mass  center  is  farther  away  from  the  elastic 
(pitch)  axis.  The  flutter  speed  increases  sharply  as  the  mass  center 

approaches  the  pitching  axis  which  is  at  J$-chord. 

4.  Flutter  Speeds  for  NACA  64A006  Airfoil 
Pitching  about  VChord  Axis  with 
Varying  ah  at  M_  = 0.7,  0.8,  and  0.85 
In  this  analysis,  the  effect  of  the  location  of  the  elastic  (pitch) 
axis  of  the  NACA  64A006  airfoil  on  flutter  speed  was  studied.  The  airfoil -air 
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REDUCED  FREQ.(kc)  FLUTTER  SPEED  (U/bw,) 


POSITION  OF  MASS  CENTERS) 

Figure  17.  Effect  of  Position  of  Mass  Center  on  Flutter  Speed  for 

NACA  64A006  Airfoil  for  Various  Frequency  Ratios  at  M = 0.8 
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INDICIAL(79x  99) 
RELAX  AT  ION(59x60) 


POSITION  OF  MASS  CENTER(xa) 

Effect  of  Position  of  Mass  Center  on  Flutter  Speed  for 
NACA  64A006  Airfoil  for  Various  Frequency  Ratios  at  M = 0.85 


Figure  18 
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mass  ratio  u was  assumed  as  100.  The  mass  center  was  assumed  to  be  constantly 
at  1/8-chord  aft  the  elastic  axis. 


Curves  for  flutter  speed  and  the  corresponding  reduced  frequency  versus 


' 

I l 


the  location  of  the  elastic  (pitch)  axis  a^  were  obtained  for  three 
different  values  of  The  elastic  axis  was  moved  from  the  20%-chord 

to  45%~chord.  These  curves  are  shown  in  Figures  19,  20,  and  21  for  Ma  = 0.7, 
0.8,  and  0.85,  respectively. 

All  three  figures  illustrate  the  same  general  trend.  The  flutter  speed 
decreases  as  both  the  elastic  axis  and  the  mass  center  move  toward  the 
trailing  edge.  The  flutter  speed  decreases  as  (o^/w^)  increases.  The 
curves  for  flutter  speed  drop  very  slightly  as  the  free  stream  Mach  number 
increases.  The  curves  for  the  reduced  frequency  demonstrate  a behavior 
opposite  to  that  of  the  curves  for  flutter  speed.  For  M^  = 0.85  and 
(o>h/t»>  ) = 0.3,  difficulties  were  encountered  in  obtaining  flutter  speed 
because  flutter  occurred  in  the  high  reduced  frequency  region. 

The  agreements  in  flutter  speed  between  the  indicial  and  relaxation 
solutions  are  good.  However,  the  flutter  speed  increases  sharply  when  the 
elastic  axis  approaches  the  leading  edge.  Thus  the  discrepancies  in  flutter 
speed  between  the  two  methods  also  become  higher  when  the  elastic  axis 
approaches  the  leading  edge. 

5^.  Effect  of  Transformation  of  Aerodynamic 

Coefficients  on  the  Flutter  Results  of 

NACA  64A006  Airfoil  at  M =0.8  and  0.85 
00  ■ — 

The  aerodynamic  coefficients  for  an  NACA  64A006  airfoil  pitching  about 
the  leading  edge  were  obtained  and  presented  in  Tables  6 and  7 for  M = 0.8 

oo 

and  0.85,  respectively.  It  is  of  interest  to  transform  these  aerodynamic 
coefficients  at  the  leading  edge  to  those  at  the  elastic  axis  (1/4-chord) 


/bu)«) 


FLUTTER  SPEED  (U/bu.) 


to  perform  a flutter  analysis.  It  is  also  of  interest  to  compare  such 

flutter  results  with  those  obtained  directly  without  such  transformation. 

Curves  for  flutter  speed  and  the  corresponding  reduced  frequency 

versus  the  airfoil-air  mass  ratio  were  obtained  for  three  different  values 

of  (w ./to  ).  The  mass  center  was  assumed  to  be  at  a distance  of  1 /8-chord 
'ho 

aft  the  elastic  (pitch)  axis  which  is  at  1/4-chord. 

The  results  obtained  on  the  basis  of  the  relaxation  method  are  shown 
in  Figures  22  and  23  for  = 0.8  and  0.85,  respectively.  In  both  figures, 
the  two  sets  of  results,  one  obtained  with  and  one  without  transformation  of 
aerodynamic  coefficients,  are  in  fairly  good  agreement. 

The  results  obtained  based  on  the  indicial  method  are  shown  in  Figures 
24  and  25  for  M = 0.8  and  0.85,  respectively.  In  these  two  figures,  the 

oo 

discrepancies  between  the  two  sets  of  results,  one  obtained  with  and  one 
without  transformation  of  aerodynamic  coefficients,  are  more  pronounced 
than  those  found  in  Figures  22  and  23. 

Such  discrepancy  may  be  attributed  to  different  numerical  errors 
introduced  in  the  two  sets  of  aerodynamic  coefficients,  one  computed  for 
pitching  axis  at  the  leading  edge  and  one  at  1 /4-chord.  It  may  also  be 
attributed  to  the  effect  of  local  angle  of  attack  which  depends  on  the 
point  of  pitching  for  the  nonlinear  case.  Such  effect  could  not  be  considered 
during  the  transformation  process. 

6.  Effect  of  Mach  Number  on  Flutter  Speed 
for  NACA  64A006  Airfoil 

Thus  far,  flutter  results  have  been  presented  only  for  given  Mach  numbers. 
It  Is  of  Interest  to  study  the  effect  of  Mach  number  on  flutter  speed  of  an 
NACA  64A006  airfoil  pitching  about  the  1 /4-chord  axis.  The  mass  center  was 
assumed  to  be  at  1/8-chord  aft  the  elastic  (pitch)  axis. 
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NACA  64A006  , 

M = 0.8 
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Figure  23.  Effect  of  Transformation  of  Aerodynamic  Coefficients  (by  Relaxa- 
tion Method)  on  Flutter  Speed  for  NACA  64A006  Airfoil  at  M = 0.85 
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Figure  24.  Effect  of  Transformation  of  Aerodynamic  Coefficients  (by  Indicia! 
Method)  on  Flutter  Speed  for  NACA  64A006  Airfoil  at  M = 0.8 
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Figure  25.  Effect  of  Transformation  of  Aerodynamic  Coefficients  (by  Indlcial 
Method)  on  Flutter  Speed  for  NACA  64A006  Airfoil  at  M = 0.85 
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Curves  for  flutter  speed  and  the  corresponding  reduced  frequency  versus 

free  stream  Mach  number  were  obtained  for  three  different  values  of  airfoil- 

air  mass  ratio  (u  * 100,  200,  and  300)  as  shown  in  Figure  26.  The  free 

stream  Mach  numbers  considered  were  0.7,  0.8,  and  0.85.  The  frequency 

ratio  (u./u)  ) was  assumed  as  0.3. 
n a 

Results  shown  were  obtained  on  the  basis  of  both  relaxation  and  indicial 
methods.  In  addition,  a curve  obtained  on  the  basis  of  linear  flat  plate 
theory  (Kernel  function  method)  for  u = 100  is  shown. 

It  is  seen  that  the  flutter  speed  reduces  as  free  stream  Mach  number 
varies  from  0.7  to  0.85.  The  rate  of  reduction  increases  as  u increases. 
Higher  M^  will  be  analyzed  based  on  the  limitions  of  the  program. 


REDUCED  FREQ. (kc)  FLUTTER  SPEED  (U/bco*) 


FREESTREAM  MACH  NUMBER  (M) 


Fig.  26.  Effect  of  Mach  number  on  Flutter  Speed  of  NACA  64A006  Airfoil  for 
Various  Frequency  Ratios. 
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SECTION  VII 


FLUTTER  RESULTS  FOR  NACA  64A010  AIRFOIL 


(a)  NACA  64A010  Airfoil  Configuration 


The  configuration  of  an  NACA  64A010  airfoil  may  also  be  described  by 


Equation  44.  There  is  no  camber  for  this  airfoil.  The  constants  for  the 


thickness  function  are:  aQ  = 0.2255645;  a^  = 0;  a2  = -0.0699163;  a3  = 
0.0664236;  a4  = -0.5415990;  and  a^  = 0.3199472.  A plot  of  the  airfoil  con- 


figuration is  given  in  Figure  27.  Again,  see  Part  IV  of  Reference  38  for 


better  coefficients. 


(b)  Aerodynamic  Coefficients 


In  this  study,  the  airfoil  again  was  assumed  to  oscillate  with  two  de- 


grees of  freedom,  plunge  and  pitch.  The  pitching  axis  was  at  the  h-ch ord. 


Thus  four  aerodynamic  coefficients  c^,  c^,  cm6,  cma  were  computed  by 


using  both  the  relaxation  method  and  the  indicial  method  for  three  different 


free  stream  Mach  numbers  0.72,  0.76,  and  0.80  and  for  various  values  of  re- 


duced frequency  k . For  the  case  M = 0.84,  however,  only  the  relaxation 

C 00 

method  was  used.  In  all  computations,  the  relaxation  method  was  based  on  a 


mesh  of  59  x 60  and  the  indicial  method  was  based  on  a mesh  of  79  x 99. 


The  results  for  the  aerodynamic  coefficients  are  shown  in  Tables  10,  11, 


12,  and  13  for  M = 0.72,  0.76,  0.80,  and  0.84,  respectively.  Between  the 


two  methods,  the  agreements  in  lift  coefficients  are,  in  general,  better  than 


the  moment  coefficients.  In  several  cases,  the  moment  coefficients  obtained 


by  the  two  methods  are  even  different  in  sign.  Also,  the  aerodynamic 


coefficients  obtained  by  the  relaxation  method  drop  off  in  high  kc  values.  This 


may  be  attributed  to  the  fact  that  the  relaxation  program  UTRANS2  has  con- 


vergence difficulty  in  higher  kc  values. 
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TABLE  10 


Aerodynamic  Coefficients  for  NACA  64A010  pitching  at  the  quarter  chord 
at  M = 0.72 


1 

1 

Reduced  Frequency 

1 

| 

k =0.05 
c 

k =( 
c 

3.10 

k =( 
c 

3.15 

k =0.20 
c 

k =0.30 
c 

■ 

Q) 

2: 

Real 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

Real 

Imag 

c. 

i 

in 

0.038 

0.452 

0.147 

0.881 

0.317 

1 .270 

0.530 

1 .600 

1 .007 

2.091 

l6 

2 

0.0 

0.0 

0.128 

0.520 

0.168 

0.780 

0.476 

1.384 

0.933 

1 .614 

Ci 

1 

9.100 

0.0 

9.030 

-0.754 

8.810 

-1 .470 

8.470 

-2.110 

8.020 

-2.650 

6.969 

-3.355 

a 

2 

9.710 

0.0 

8.933 

-1.171 

8.197 

-1.851 

6.679 

-2.289 

5.376 

-1 .993 

c 

1 

0.0 

0.0 

0.000 

-0.014 

-0.001 

-0.028 

-0.001 

-0.042 

-0.002 

-0.056 

-0.005 

-0.084 

m6 

2 

0.0 

0.0 

-0.001 

-0.015 

0.007 

-0.021 

0.029 

-0.051 

0.038 

-0.091 

Cm 

1 

-0.282 

0.0 

-0.282 

0.003 

-0.282 

0.005 

-0.281 

0.008 

-0.281 

0.010 

-0.281 

0.015 

Tna 

2 

-0.240 

0.0 

-0.231 

-0.062 

-0.237 



-0.122 

■ 

-0.253 

-0.228 

-0.265 

-0.310 

Method  1:  Indicial  Method  (79  x 99) 
Method  2:  Relaxation  Method  (59  x 60) 
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TABLE  11 


Aerodynamic  Coefficients  for  NACA  64A010  pitching  at  the  quarter  chord  at 


Reduced  Frequency 


kc=°.° 


k =0.05 
c 


k =0.10 
c 


k =0.15 
c 


kc=0.20 


k =0.25 
c 


Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

0.068 

0.519 

0.252 

0.976 

0.503 

1 .330 

0.762 

1 .590 

0.985 

1.770 

0.106 

0.495 

0.273 

0.860 

0.501 

1 .100 

0.668 

1.258 

0.850 

1.370  ! 

10.380 

-1 .360 

9.760 

-2.520 

8.890 

-3.350 

7.930 

-3.810 

7.070 

-3.940 

9.766 

-1 .682 

8.650 

-2.485 

7.620 

-2.769 

6.724 

-2.711 

5.925 

-2.553 

0.000 

-0.019 

-0.001 

-0.037 

-0.003 

-0.056 

-0.005 

-0.075 

-0.008 

-0.093 

0.001 

-0.014 

0.007 

-0.027 

0.013 

-0.043 

0.020 

-0.058 

0.025 

-0.079 

-0.374 

0.006 

-0.374 

0.013 

-0.373 

0.019 

-0.373 

0.025 

-0.372 

0.032 

-0.264 

-0.068 

-0.271 

-0.134 

-0.283 

-0.194 

-0.293 

-0.249 

-0.305 

-0.290 

Method  1:  Indiclal  Method  (79  x 99) 
Method  2:  Relaxation  Method  (59  x 60) 


TABLE  12 


Aerodynamic  Coefficients  for  NACA  64A010  pitching  at  the  quarter 
chord  at  M = 0.80 


Reduced  Frequency 


k =0.05 
c 

k =0.10 
c 

k =0.15 
c 

k =0.20 
c 

Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

0.112 

0.638 

0.408 

1.160 

0.797 

1 .500 

1 .160 

1 .670 

0.143 

0.557 

0.394 

0.918 

0.639 

1.117 

0.860 

1 .235 

12.750 

-2.230 

11 .610 

-4.080 

10.030 

-5.310 

8.350 

-5.820 

11.109 

-2.673 

9.211 

-3.621 

7.702 

-3.647 

6.510 

-3.352 

0.000 

-0.028 

-0.001 

-0.057 

-0.002 

-0.086 

-0.003 

-0.116 

0.001 

-0.021 

0.005 

-0.041 

0.009 

-0.063 

0.012 

-0.088 

-0.564 

0.005 

-0.567 

0.009 

-0.572 

0.013 

-0.582 

0.017 

-0.409 

-0.060 

-0.407 

-0.126 

-0.411 

-0.192 

-0.417 

-0.244 

TABLE  13 
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Aerodynamic  Coefficients  for  NACA  64A010  pitching  at  the  quarter 
chord  at  M = 0.84 


1 

Reduced  Frequency 

k =0.0 
c 

k =0.05 
c 

k =0.10 
c 

k =0.15 
c 

k =0.20 
c 

Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

Real 

Imaq 

c 

*6 

mm 

0.126 

0.469 

0.324 

0.791 

0.530 

0.945 

0.695 

1.017 

10.914 

0.0 

9.256 

-2.177 

7.746 

-2.893 

6.466 

-2.993 

5.267 

-2.771 

Cm6 

B 

0.0 

-0.019 

-0.075 

-0.049 

-0.130 

-0.087 

-0.159 

-0.126 

-0.176 

Cm 

a 

-1 .719 

-1.474 

0.283 

-1 .266 

0.364 

-1 .094 

0.370 

-0.932 

0.363 

Relaxation  Method  (59  x 60) 


The  distributions  of  the  steady  pressure  coefficient  along  the  airfoil 
obtained  by  using  both  STRANS2  and  LTRAN2  for  the  NACA  64A010  airfoil  at 
= 0.72,  0.76,  0.80,  and  0.84  are  shown  in  Figures  28,  29,  30,  and  31, 
respectively.  It  should  be  noted  that  steady  solutions  in  both  STRANS2 
and  LTRAN2  are  based  on  the  relaxation  method.  The  discrepancies  between  the 
two  solutions  may  lie  in  the  difference  in  mesh  and  the  difference  in  the 
procedure  of  treating  the  shock. 

(c)  Flutter  Results 

In  the  flutter  analysis  of  the  NACA  64A010  airfoil,  the  standard  U-g 
method  was  used.  The  effects  of  four  parameters  were  considered:  airfoil- 
air  mass  ratio  u;  position  of  the  mass  center  xq;  position  of  the  elastic 
axis  a^;  and  the  plunge  - pitch  frequency  ratio  (<^/wa).  In  all  the  following 
cases  the  values  for  the  radius  of  gyration  r^  and  the  reference  frequency 
u>r  were  assumed  as  0.5  and  1.0,  respectively. 

1 . Flutter  Speeds  for  NACA  64AQ10  Airfoil  Pitching  about  fr-Chord  Axis  with 
Varying  u at  M = 0.72,  0.76,  and  0.8 

The  effect  of  airfoil  - air  mass  ratio  was  first  studied.  The  NACA 

64A010  airfoil  was  assumed  to  pitch  about  the  elastic  axis  located  at  \- 

chord  (ah  =-0.5).  The  mass  center  was  assumed  to  be  at  1/8-chord  aft  the 

elastic  axis  (xa  = >4)  or  at  37.5%-chord  from  the  leading  edge.  The  range 

of  u considered  was  between  25  to  300.  The  frequency  ratios  (u>.  / u ) con- 

n a 

sldered  were  0.1,  0.2,  0.3,  and  0.4 

Figure  32  shows  a plot  of  the  nondimensional  flutter  speed  (U/bu  ) and 

a 

the  corresponding  reduced  frequency  ( kc ) versus  airfoil-air  mass  ratio  at 
= 0.72  and  = 0.1,  0.2,  0.3,  and  0.4,  respectively.  It  is  seen  that 

at  <*>h/«a  = 0.1,  the  flutter  speed  Increases  as  v increases.  Such  increasing 
rate  reduces  as  the  frequency  ratio  a Increases.  The  curves 
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Figure  28.  Distribution  of  Steady  Pressure  Coefficients  for 
NACA  64A010  Airfoil  at  M = 0.72 
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Figure  31.  Distribution  of  Steady  Pressure  Coefficients  for 
NACA  64A010  Airfoil  at  M = 0.84 
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Figure  32.  Effect  of  Mass  Ratio  on  Flutter  Speed  for  NACA  64A010  Airfoil 
for  Various  Frequency  Ratios  at  M = 0.72 


for  flutter  speed  are  lower  for  higher  frequency  ratios 

The  behavior  of  the  curves  for  reduced  frequency  Is  opposite  to  that  of 
the  curves  for  flutter  speed.  In  the  lower  y range,  the  reduced  frequencies 
that  cause  flutter  are  higher  than  the  valid  limit  of  the  present  compu- 
tational methods. 

Results  similar  to  those  given  in  Figure  32  are  presented  In  Figures  33 
and  34  for  = 0.76  and  0.8,  respectively.  Similar  trends  are  observed.  The 
curves  for  flutter  speed  drop  slightly  as  the  free  stream  Mach  number  Increases 
It  is  interesting  to  see  that  for  higher  free  stream  Mach  number  and  higher 
frequency  ratio,  the  variation  of  y value  has  less  effect  on  the  flutter 
speed. 

It  may  be  pointed  out  that  for  all  the  three  cases  (Mro  = 0.72,  0.76,  and 

0. 8),  the  flutter  results  obtained  by  using  both  the  relaxation  method  and 
the  indicial  method  are  in  fairly  good  agreement. 

1.  Flutter  Speeds  for  NACA  64A010  Airfoil  Pitchinq  about  Jj-Chord  Axis  with 

Varying  x at  M = 0.72,  0.76,  and  0.8 

In  this  case,  the  effect  of  the  position  of  the  mass  center,  defined  as 
a distance  xab  aft  the  elastic  axis,  on  the  flutter  speed  and  the  corres- 
ponding reduced  frequency  is  studied. 

Again,  the  NACA  64A010  airfoil  was  assumed  to  pitch  about  the  ?*-chord 

axis  (a^  = -0.5).  The  airfoil  - air  mass  ratio  y was  assumed  as  100.  The 

frequency  ratios  (u>^/u)a)  considered  were  0.1,  0.2,  0.3,  and  0.4.  The  mass 

center  was  assumed  to  move  from  the  elastic  axis  at  %-chord  to  %-chord 

(from  x = 0 to  x =0.5). 
a a 

Plots  of  the  nondimensional  flutter  speed  and  the  corresponding  reduced 
frequency  versus  the  position  of  mass  center  x_  are  shown  in  Figures  35,  36, 
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Figure  34.  Effect  of  Mass  Ratio  on  Flutter  Speed  for  NACA  64A010 
Airfoil  for  Various  Frequency  Ratios  at  M = 0.8 
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Figure  35.  Effect  of  Position  of  Hass  Center  on  Flutter  Speed  for  NACA 
64A010  Airfoil  for  Various  Frequency  Ratios  at  M * 0.72 
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Figure  37.  Effect  of  Position  of  Mass  Center  on  Flutter  Speed  for 

NACA  64A010  Airfoil  for  Various  Frequency  Ratios  at  M * 0.8 


and  37  for  free  stream  Mach  numbers  of  0.72,  0.76,  and  0.8,  respectively. 

All  three  figures  Illustrate  the  same  general  trend.  The  flutter  speed 
Increases  steadily  as  the  mass  center  moves  away  from  the  mid-chord  toward 
the  elastic  axis.  However,  the  flutter  speed  increases  sharply  as  the  mass 
center  becomes  closer  to  the  elastic  axis.  The  flutter  speed  decreases  as 
the  frequency  ratio  increases.  The  curves  for  flutter  speed  drop  very 
slightly  as  the  free  stream  Mach  number  increases. 

The  curves  for  the  reduced  frequencies  kc  assume  a general  trend  oppo- 
site to  that  of  the  curves  for  the  flutter  speed. 

The  agreement  between  the  two  sets  of  curves,  one  obtained  based  on 
the  relaxation  method  and  the  other  on  the  Indicial  method,  is  quite  good 
for  all  the  three  free  stream  Mach  numbers.  It  may  be  pointed  out,  however, 
such  agreement  becomes  poor  when  the  flutter  speed  increases  sharply  at  the 
lower  xa  range. 

.3.  Flutter  Speeds  forNACA  64A010  Airfoil  Pitching  about  fr-Chord  Axis  with 
Varying  a^  at  Mx  = 0.72,  0,76,  and  0.8 


In  this  case,  the  effect  of  the  location  of  the  elastic  (pitch) 
axis  of  the  NACA  64A010  airfoil  on  the  flutter  speed  and  the  corresponding 
reduced  frequency  was  studied.  The  airfoil  - air  mass  ratio  y was  assumed 
as  100.  The  mass  center  was  assumed  to  be  always  1 /8-chord  aft  the 
elastic  axis.  The  frequency  ratios  considered  were  0.1,  0.2,  0.3,  and  0.4. 

The  elastic  axis  was  assumed  to  move  from  20%-chord  to  45%-chord  (from  ah  = 

-0.6  to  ah  * -0.1). 

Plots  of  the  nondimenslonal  flutter  speed  and  the  corresponding  reduced 
frequency  versus  the  position  of  the  elastic  axis  a^  are  shown  in  Figures 
38,  39,  and  40  for  free  stream  Mach  numbers  of  0.72,  0.76,  and  0.8,  respectively. 
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All  three  figures  illustrate  the  same  general  trend.  The  flutter  speed 
decreases  as  both  the  elastic  axis  and  the  mass  center  move  toward  the 
trailing  edge.  The  curves  for  flutter  speed  drop  as  the  frequency  rat.o 
(»h/»0)  increases  =»nd  they  also  drop  slightly  as  the  free  stream  Mach 
number  increases  from  0.72  to  0.80.  The  curves  for  the  reduced  frequency 
present  a general  trend  opposite  to  that  of  the  curves  for  the  flutter 
speed.  At  M = 0.76  and  0.8,  it  was  not  possible  to  obtain  flutter  speed 

oo 

for  higher  frequency  ratio  such  as  = 0,4  because  the  flutter  occurred 

in  the  high  reduced  frequency  region. 

The  flutter  speed  curves  obtained  by  using  both  the  indicial  and  the 
relaxation  method  are  in  good  agreement  in  all  the  three  figures.  However, 
due  to  the  fact  that  the  flutter  speed  increases  sharply  when  the  elastic 
axis  approaches  the  leading  edge,  the  discrepancies  in  flutter  speed  between 
the  two  methods  also  become  higher  when  the  elastic  axis  approaches  the  lead 
ing  edge. 


4.  Effect  of  Mach  Number  on  Flutter  Speed  for  NACA  64A010  Airfoil 


In  this  section,  a set  of  nine  curves  has  been  presented  to  show  the 

effects  of  airfoil  - air  mass  ratio  y,  position  of  mass  center  x , position 

a 

of  elastic  axis  a^,  and  frequency  ratio  (w^/w^)  on  the  flutter  speed  of  a 
NACA  64A010  airfoil  at  three  different  free  stream  Mach  numbers  (0.72,  0.76, 
and  0.8).  It  is  of  interest  to  study  the  effect  of  the  free  stream  Mach 
number  on  the  flutter  speed  of  such  an  airfoil  pitching  about  the  h- chord  axis 
In  this  study,  the  mass  center  was  assumed  to  be  at  a distance  of  1 /8-chord 
aft  the  elastic  (pitch)  axis  (xa  = 0.25).  The  frequency  ratio  (u^/u^) 
was  assumed  as  0.3. 
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The  results  were  obtained  as  plots  of  flutter  speed  and  the  corres- 
ponding reduced  frequency  versus  free  stream  Mach  number  for  three  diff- 
erent values  of  airfoil  - air  mass  ratio  (y  s 100,  200,  and  300)  as  shown 
In  Figure  41.  The  free  stream  Mach  numbers  considered  were  0.72,  0.76,  and 


The  flutter  speed  decreases  as  the  free  stream  Mach  number  Increases 
from  0.72  to  0.8.  The  rate  of  decrease,  however,  decreases  as  the  airfoil 


air  mass  ratio  y becomes  smaller. 

The  two  sets  of  curves,  one  obtained  by  using  the  relaxation  method  and 
the  other  by  the  indiclal  method,  show  the  same  general  trend. 


REDUCED  FREQ.(kc)  FLUTTER  SPEED 


SECTION  VIII 
CONCLUDING  REMARKS 

The  purpose  of  this  research  was  to  study  the  flutter  behavior  of  two- 
dimensional  and  two  degree-of- freedom  airfoils  In  small  disturbance,  un- 
steady transonic  flow.  Two  computer  programs,  STRANS2  and  UTRANS2  based 
on  the  harmonic  or  relaxation  method  and  LTRAN2  based  on  the  Indlclal  and 
time- Integration  methods,  were  used  to  compute  the  aerodynamic  coefficients 
for  flutter  analysis.  An  NACA  64A006  and  an  NACA  64A010  airfoil  were  selected 
for  this  study.  A flat  plate  was  also  selected  In  the  study  In  order  to 
compare  the  flutter  solutions  based  on  the  linear  flat  plate  theory  with 
those  based  on  the  present  two  computer  programs  at  a low  free  stream  Mach 
number  of  0.7. 

In  the  flutter  analysis,  the  effects  of  four  parameters  were  considered: 

(1)  the  airfoil  - air  mass  ratio  p;  (2)  the  position  of  the  mass  center  xq; 

(3)  the  position  of  the  elastic  (pitch)  axis  ah;  and  (4)  the  plunge  to 
pitch  frequency  ratio  (u>h/u)a).  The  Mach  numbers  considered  were  0.7,  0.8, 
and  0.85  for  the  NACA  64A006  airfoil  and  0.72,  0.76,  and  0.8  for  the  NACA 
64A010  airfoil. 

As  a result  of  this  study,  the  following  concluding  remarks  may  be  made: 

(1)  Although  the  time- Integration  method  Is  more  accurate  than  the 
relaxation  method.  It  requires  relatively  more  computing  time.  Thus  the 
Indlclal  method  was  used  to  compute  all  the  present  aerodynamic  coefficients. 

In  order  to  compute  all  the  aerodynamic  coefficients  at  M^  * 0.76  for  the 
NACA  64A010  airfoil  for  five  different  values  of  k as  shown  In  Table  11, 

C 

the  Indlclal  method  required  approximately  2.1  hours  central  processing  (CP)  time 
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on  a CDC  6500  computer  whereas  the  relaxation  method  required  4.7  hours. 


If  the  time- integration  method  Is  used  for  the  same  mesh  as  the  indlclal 
method,  approximately  7 hours  CP  time  will  be  needed.  In  the  present  com- 


putations, the  finite  difference  mesh  plays  a very  important  role. 


(2)  At  a low  free  stream  Mach  number  of  0.7,  the  results  for  the 


steady  pressure  curves,  the  unsteady  pressure  curves,  the  aerodynamic  co- 


efficients, and  the  flutter  speeds  obtained  by  using  the  computer  programs 


STRANS2,  UTRANS2,  LTRAN2,  and  the  linear  flat  plate  theory  program  based 


on  the  Kernel  function  method  agree  reasonably  well. 


(3)  The  steady  solutions  of  both  programs  STRANS2  and  LTRAN2  are 


based  on  the  same  relaxation  method.  Thus  the  steady  pressure  distribution 


curves  obtained  by  using  both  programs  are  in  good  agreement  in  the  low 


Mach  numbers  that  result  in  no  shocks  and  present  some  difference  in  the 


higher  Mach  numbers  that  result  in  shocks.  This  difference  may  be  due  to 


the  difference  in  the- mesh  size  and  in  the  manner  the  shocks  are  treated 


in  the  two  programs.  Because  in  both  the  relaxation  and  indicial  methods 


the  unsteady  computation  depends  on  the  steady  solution  as  input,  the 


differences  in  steady  results  contribute  to  the  disagreement  in  the  unsteady 


results. 


(4)  The  aerodynamic  coefficients  obtained  by  using  both  the  relaxation 


and  the  indicial  methods  are  compared  in  Tables  1-4,  6-7,  and  10-12.  The 


lift  coefficients  agree  better  than  the  moment  coefficients.  This  is  due 


to  the  fact  that  the  lift  coefficients  were  obtained  by  merely  integrating 


the  area  under  the  pressure  curves  whereas  the  moment  coefficients  were  ob- 


tained by  integrating  the  area-moment  urder  the  pressure  curves.  The  former 


integration  resulted  in  less  error  than  the  latter.  Furthermore,  the  mag- 


nitudes of  the  moment  coefficients  depend  on  the  location  of  the  pitching 
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axis.  They  are  relatively  small  when  the  pitching  axis  Is  at  the  %-chord. 

In  that  case,  the  percentage  differences  In  moment  coefficients  between  the 
two  methods  could  be  magnified  and  may  not  be  significant. 

(5)  Due  to  the  low  frequency  approximation  used  In  LTRAN2,  the  re- 
duced frequencies  k£  considered  In  most  of  the  cases  were  not  higher  than 
approximately  0.2.  Although  the  low  frequency  approximation  was  not  used 

In  UTRANS2,  convergence  could  not  be  obtained  at  higher  reduced  frequencies. 
In  this  study,  the  reduced  frequencies  considered  were  never  higher  than 
0.3. 

(6)  The  curves  for  the  flutter  speed  and  the  corresponding  reduced 
frequency  obtained  by  using  both  the  relaxation  and  the  Indlclal  methods 
are.  In  general.  In  fairly  good  agreement. 

(7)  The  main  difficulty  In  obtaining  flutter  results  was  In  the 
upper  limitation  of  the  reduced  frequency  k . In  all  the  cases,  the  search 
for  flutter  points  was  below  the  k£  values  of  0.25  for  LTRAN2  and  0.3  for 
UTRANS2.  Due  to  such  limitations,  It  was  not  possible  to  obtain  flutter 
results  for  the  airfoil  - air  mass  ratios  p less  than  approximately  25  and 
for  the  frequency  ratios  <i>h/wa  higher  than  approximately  0.4. 

(8)  The  general  trend  of  the  present  curves  for  flutter  speed  are 
similar  to  the  curves  In  subsonic  region. 

(9)  The  effect  of  each  of  the  parameters  p,  xa,  ah,  and  M* 

Is  presented  and  discussed  In  detail  In  Sections  VI  and  VII. 

(10)  Based  on  the  limited  Mach  numbers  considered.  It  was  observed 
that  the  flutter  speed  reduces  as  the  Mach  number  Increases.  But  the  rate 
of  reduction  decreases  as  the  airfoil  - air  mass  ratio  decreases.  More 
general  conclusions  may  be  made  If  more  and  higher  Mach  numbers  were  con- 
sidered. However,  the  present  computer  programs  are  not  valid  for  higher 
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transonic  Mach  numbers  because  the  Invlscld  flow  theory  Is  assumed. 

(11)  This  analysis  may  provide  a systematic  procedure  for  the  transonic 
flutter  analysis  of  airfoils.  The  study  may  further  provide  a review  of  the 
limitations  and  capabilities  of  the  two  existing  transonic  codes.  The 
present  results  may  give  some  Insight  Into  the  transonic  flutter  behaviors 
of  the  airfoils.  The  present  results  may  also  be  used  as  a comparative  basis 
by  the  experimentalists  and  other  analysts  who  work  In  the  same  field. 

(12)  In  order  to  make  more  rigorous  conclusions,  the  study  of  more 
airfoils,  especially  the  supercritical  ones,  is  needed. 

(13)  For  a more  complete  understanding  of  the  transonic  flutter  behaviors 
of  the  airfoils  In  wider  ranges  of  parameters,  the  present  limitations  on 
very  low  reduced  frequency  and  low  transonic  Mach  numbers  must  be  removed. 
Improvement  In  transonic  codes  to  allow  for  higher  values  of  k£  and  Is 
needed . 

(14)  To  be  more  realistic,  the  structural  flexibility  of  the  airfoil 
should  be  considered.  This  can  be  accounted  for  by  using  curved  beam-type 
finite  elements. 

(15)  In  order  to  study  more  practical  cases  of  full  wing  flutter,  the 
extension  of  the  present  2-D  transonic  code  to  account  for  the  3-D  case  Is 
needed . 
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